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In this work we exploit the integrability of the two-lead Anderson model to compute 
transport properties of a quantum dot, in and out of equilibrium. Our method combines the 
properties of integrable scattering together with a Landauer-Buttiker formalism. Although 
we use integrability, the nature of the problem is such that our results are not generically 
exact, but must only be considered as excellent approximations which nonetheless are valid 
all the way through crossover regimes. 

The key to our approach is to identify the excitations that correspond to scattering 
states and then to compute their associated scattering amplitudes. We are able to do so 
both in and out of equilibrium. In equilibrium and at zero temperature, we reproduce 
the Friedel sum rule for an arbitrary magnetic field. From this we compute exactly the 
behaviour of the zero temperature linear response conductance as a function of both the 
gate voltage and the field. We also study transport quantities requiring knowledge of 



scattering states away from the Fermi surface. We compute the linear response conductance 
at finite temperature at the symmetric point of the Anderson model, and reproduce Costi 
et al.'s numerical renormalization group computation of this quantity. We then explore 
the out-of-equilibrium conductance for a near-symmetric Anderson model, and arrive at 
quantitative expressions for the differential conductance, both in and out of a magnetic 
field. We reproduce the expected splitting of the differential conductance peak into two in 
a finite magnetic field, H. We determine the width, height, and position of these peaks. 
In particular we find for H ^> T^, the Kondo temperature, the differential conductance 
has maxima of e 2 /h occuring for a bias V close to but smaller than H. The nature of 
our construction of scattering states suggests that our results for the differential magneto- 
conductance are not merely approximate but become exact in the large field limit. 
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1. Introduction 

The Kondo effect is a cynosure of modern condensed matter physics. Due to the 
strongly coupled nature of its IR fixed point, understanding its low energy behaviour has 
proven a major theoretical challenge. Typically the phenomena refers to the interaction of 
isolated magnetic impurities in a bulk metal. However in the last several years the exper- 
imental study of single magnetic impurities has moved to a new arena, that of quantum 
dots connected to external leads M H] B [Q 0. In analogy to the traditional realization 
of the Kondo effect, the leads serve as the bulk metal and the dot as the magnetic impurity. 
The appearance of the Kondo effect in this new setting has brought a new set of theo- 
retical challenges: how to compute transport quantities that form the main experimental 
signatures of these systems. 

Quantum dots come in at least two forms. Semi-conductor quantum dots M O M 0] are 
a product of the continuing project of the miniaturization of solid state devices. They 
are fabricated by confining electrons in a two dimensional electron gas (2DEG) within a 
GaAs/AlGaAs structure using a combination of metallic gates. The region to which the 
electrons are confined is small enough that its energy levels may be considered discrete. 
The dot is connected to source and drain contacts (the two leads). Schematically the 
quantum dot can be pictured as follows: 






Figure 1: A schematic of the quantum dot. 
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The source and drain can be held at any relative voltage, thus enabling the study of both 
linear response and out-of-equilibrium transport quantities. Beyond the gates that serve to 
confine the electrons of the 2DEG, additional gates can be deposited on the GaAs/AlGaAs 
heterostructure. Such gates capacitively couple to the quantum dot through a gate voltage, 
V g , thus allowing the chemical potential of the dot to be adjusted. This has two important 
consequences. By adjusting the gate voltage, one can tune the number of electrons in the 
confined region to be odd so that there is a single unpaired electron. The presence of an 
unpaired electron allows the appearance of Kondo-like physics. Moreover by tuning V g , 
the unpaired electron's chemical potential can be adjusted thus controlling the scale, Tk, 
where Kondo physics sets in. Such systems are thus said to possess a "tunable" Kondo 
effect @. 

Quantum dots may also be fabricated from metallic carbon nanotubes || . By deposit- 
ing metallic leads on top of a small section of a carbon nanotube, an effective quantum 
dot is made. Like their semiconductor counterparts, these dots are tunable: gates may be 
added to the semiconductor substrate upon which the nanotube/leads lay. Semiconductors 
dots typically carry 10-100 electrons. Nanotube dots, in contrast, have many thousands of 
electrons and yet still exhibit Kondo-like physics. 

The transport quantities that lie at the focus of the experimental study of a tunable 
Kondo effect in quantum dots have been measured under a variety of conditions. Con- 
ductances of the lead-dot system have been determined both in and out of equilibrium 
and both in and out of the presence of a magnetic field. Remarkably this wide variety of 
experimental phenomena is thought to be described by a conceptually simple theory: the 
Anderson model. 

The Anderson model is fashioned from a chain of non-interacting spinful fermions (the 
leads) connected via hopping to a single site impurity on which alone Coulomb repulsion 
is present. On a lattice, the Hamiltonian reads (with no bias and H=0): 

H = -t(c\ ia Ci +1)0t + h.c.) + Un^ridi 

(1.1) 

+ {Vi( c -i,a d a + h.c) + V 2 (c\ a d a + h.c)) + e d (n d ^ + n d [). 

a. 

Here c^/c a are the lead electron creation/destruction operators, d' /d the dot electron 
creation/destruction operators, and n d = d^d. The dot is considered to reside at x = 0. 
The index a indicates spin species. The interaction on the dot is present in the term 
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Un^Tidi- Although deceptively simple, the presence of a non-zero U makes the problem 
many-body with all of its manifold complications. 

U is pictured in Figure 1 and represents the charging energy incurred when an electron 
is added to the dot. Roughly it can be estimated as 

u = w + Ae ' (L2) 

where C is the capacitance of the dot and Ae is the dot's energy level spacing. For the 
experiments at hand, U ~ lmeV. The counterpart of the gate voltage, V g , in the above 
Hamiltonian is Ed, the parameter that controls the chemical potential of the electrons on 
the dot. By adjusting Ed, the number of electrons on the model dot can be varied from zero 
to two. (Although there may be a large number of electrons on the actual dot, the concern 
here both theoretically and experimentally is of electrons in the highest occupied energy 
level.) The final pair of parameters, Vi^, measure the height of the tunnel barriers between 
the leads and the dot. In general, they differ between the two leads. The quantity, T, (see 
Figure 1) measuring the width of the level resonance that arises through the interaction of 
the leads and the dot is given in terms of the tunnel barrier heights to be Y = (V^ 2 + V 2 2 )/2. 
Typically T/U -C 1 in the experiments of concern [0@@@. 

Although the above described experimental results have come relatively recently, the 
theoretical study of transport through impurities is much older. Appelbaum and Anderson 
both studied conductance anomalies present in tunnel junctions due to the presence of 
magnetic impurities 0. However their efforts were perturbative in nature and did not 
describe the strong coupling nature of the Kondo effect. More recently, Ng and Lee 0, 
studied the linear response conductance both in and out of a magnetic field using the 
Friedel sum rule. The Friedel sum rule relates the scattering phase of the electrons at the 
Fermi surface to the average number of electrons sitting on the dot. However the Friedel 
sum rule is useful only in determining the linear response conductance, a consequence of 
the rapid variation in the scattering phase as one moves away from the Fermi surface. 
In contrast to the linear response conductance where the Friedel sum rule is an exact 
relationship, the techniques used to determine out-of-equilibrium transport are limiting in 
nature. In one approach, a non-crossing approximation (NCA) || |J was employed. The 
NCA approach has drawbacks. In order to implement the associated use of slave bosons, 
one must take U = oo. Moreover, NCA is in some sense a large N approximation where N 
is the number of spin degrees of freedom of the impurity (in this case N=2). It is known to 
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be remarkably accurate in computing thermodynamics. However it is less accurate when 
it comes to transport quantities (~ 15% errors @) due to its less accurate prediction of 
behaviour right at the Fermi surface. And as such these difficulties render it unusable in 
non-zero magnetic fields ||. In another approach to computing non-equilibrium properties, 
a clever combination of the analysis of the equations of motion with perturbation theory 
was employed to study the differential magneto-conductance [|10| . However the truncation 
of the equations of motion necessary to perform the analysis in this work is in some 
sense an uncontrolled approximation. The authors of indicate that their methodology 
underestimates the magnitude of the differential conductance. Another set of approaches 
have relied upon perturbation theory |j| [|13| . As with the results of |§ , perturbation 
theory requires relatively small U (Coulomb repulsion) or alternatively, temperatures far 
in excess of the Kondo temperature, and so presumably can access, at best, qualitative, 
not quantitative, features of the strongly coupled physics found in the Kondo regime of 
quantum dots at low temperatures. 

These inherent difficulties with the out-equilibrium Anderson were circumvented in 



the study of a non-equilibrium Kondo impurity at its Toulouse point fll4}] . At this point, 
the model can be mapped to a system of non-interacting fermions, thus permitting an ex- 
act solution. It is unclear, however, how the Toulouse limit affects the underlying physics. 
Although the ordinary Kondo model shares the same IR fixed point as its Toulouse coun- 
terpart, we are interested in part in physics for large applied field, bias, and temperature, 
that is, in physics far away from this fixed point. 

Given the limitations of these methods, one cannot help but notice that the Anderson 
model is exactly solvable. Indeed this integrability has already been exploited through 



Bethe ansatz solutions to compute thermodynamic quantities [|T6[ such as the specific 
heat and magnetic susceptibility. But what of transport quantities? A limited attempt 
to deduce information about transport properties from the Bethe ansatz solution of the 



Kondo model was made recently ||17|| . There the equilibrium impurity density of states that 
arises from the Bethe ansatz was studied. In general, the impurity density of states coming 
from the Bethe ansatz is unrelated to the spectral density of states arising from the dot 
correlator, Im(dS). Indeed at zero temperature and zero field, it is clear the two are much 
different quantities (the heights of zero energy peaks in both quantities are controlled 
by far different energy scales). But in methodology of |T7[], it is this latter quantity, 
Im(dd^), that is directly related to transport |T^]0[^O|. Moreover, the context of their 
computation, as determined in fT2|| [P| [p!U[ , demands that the non- equilibrium properties of 
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Im(dd^) be computed. Given the general unavailability from integrability of information 
about correlators such as lm{dcft), a different approach is needed to extract transport 
properties from the exact solvability of the model. Here we advocate a Landauer-Biittiker 
approach to transport, and so are instead faced with the task of identifying scattering 
states in the context of integrability. 

The key feature of an integrable system is the exact knowledge of eigenfunctions of 
the fully interacting Hamiltonian. In turn there is a well-defined notion of elementary 
excitations. In particular these excitations have an infinite lifetime: integrability forbids 
any decay processes from occurring. This arises from the infinite series of non-trivial 
conservation laws in the model. In some sense an integrable system is a superior version 
of a Fermi liquid. 

In the Anderson model, there is such a set of excitations, as detailed in Sections 2 
and 4. They are not on the face of it, however, particularly electronic. And if we are to 
understand the transport of the sea of electrons in the attached leads, we necessarily need 
scattering states which carry the quantum numbers of an electron. Rather the excitations 
divide into separate spin and charge sectors. The closest they come to being electronic is 
in bound states between excitations which can be thought of as bound states of electrons. 
This is not so unnatural. If one were to bosonize the Anderson model, one would find that 
the degrees of freedom separate into spin and charge bosons. But this is only one problem 
with the excitations arising from integrability. These excitations, as explained in Section 
2, are a combination of degrees of freedom in both of the leads connected to the dot. And 
it is the case that this entanglement cannot always be simply reversed. 

And so there is the difficulty. The scattering states are not necessarily electronic in 
nature and not confined to a single lead. Only if one can understand electronic excitations 
in an individual lead can one hope to make sense of scattering amplitudes off the dot. It 
is these two facts that have prevented the integrability of the Anderson model from being 
applied to transport quantities up to now. 

However we have managed to circumvent these problems in a number of cases. In 
particular we have successfully described both scattering states at the Fermi surface for 
generic values of U, Ed, and T, and scattering at finite energies at the symmetric point of the 
Anderson model, U = —Ed/2. There we argue that by correctly gluing together a spin and 
charge excitation, we are able to form an excitation that is electronic in nature. Moreover 
the excitations are such that one can understand them in terms of the individual leads and 
so compute reflection and transmission amplitudes of the excitation off the dot. We do so 
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in an argument akin to that used by N. Andrei [18] in computing the magnetoresistance 
in the Kondo model. There he argues that the scattering phase of an excitation can be 
identified with its impurity momentum. In turn this momentum is related to the impurity 
density of states as it appears in the Bethe ansatz and so can be directly computed. 

We now turn to how we use these excitations to compute transport quantities. All such 
quantities could be expressed in terms of the scattering of asymptotically free electrons (i.e. 
electrons in the attached leads) off the quantum dot. However such scattering, in general, 
is not particular simple. In general away from the Fermi surface such scattering will be 
inelastic and involve particle-hole production. We, however, recast the density of states 
of asymptotically free fermions in terms of the integrable excitations we have identified. 
Because of their integrability, their scattering is simple: their character does not change in 
scattering through the dot and their transport can be described individually: they scatter 
one-by-one through the dot. 

It is however unlikely that the integrable electronic excitations we use in computing 
transport properties provide exact results in all cases - the issues involved here are subtle 
and will be discussed in detail in the next sections. In particular, it is unlikely the high 
energy limit of the excitations we construct are entirely confined to a single lead. However 
this methodology successfully passes a number of tests. The first test of our method comes 
in proving the Friedel sum rule. The Friedel sum rule relates the occupancy of spin t / J. 
electrons on the quantum dot, n^/i, to the scattering phase of an electron of the same 
spin, S e yi, at the Fermi surface: 

S eVl = 27Tn Vl- (1-3) 

It thus relates a dynamic quantity to a thermodynamic quantity. With this in hand, 
previous works have computed the linear response conductance from the knowledge of this 



occupancy, at least at H = [[19] [pD[ . However such works do not make any attempt to 
explicitly identify the excitations that scatter according to the Friedel sum rule. Here we 
do so. We show that the scattering phase of the excitations we have identified to be the 
same as that predicted by the Friedel sum rule both in and out of a magnetic field. As we 
have reproduced the Friedel sum rule, we can say that the excitations we have identified 
coincide exactly with the free fermions at the Fermi surface. 

Now while the Friedel sum rule only deals with excitations at the Fermi surface, our 
methods goes beyond excitations directly at the Fermi surface, at least near the symmetric 
point of the Anderson model. To determine whether the excitations we have identified 



7 



together with their associated scattering amplitudes provide a complete solution of the 
problem, we compute the linear response conductance at finite temperature. At finite 
temperature, the excitations needed for the linear response conductance using a Landauer- 
Biittiker formalism exist over a range of energies. In particular, we compute the linear 
response conductance at the symmetric point (-U/2 = ea) of the Anderson model as a 
function of temperature, T, and compare it to Costi et al.'s f21j numerical renormalization 
group (NRG) computation of this quantity and find excellent agreement. We thus are able 
to conclude that using our excitations away from the Fermi surface is a valid procedure. 

Although our finite temperature computation suggests we have correctly identified 
the low (but finite) energy excitations at the symmetric point, we do not claim that 
our result is exact. Again our inability to make this claim hinges on the question of 
whether the integrable excitations we construct are entirely confined to a single lead and so 
make appropriate scattering states. Moreover we know that our prescription for scattering 
fails once we leave the Kondo regime where approximately one electron sits on the dot 
and enter the mixed valence regime of the Anderson model. This again suggests that 
near the symmetric point, our methods are merely highly accurate. The situation here is 
not dissimilar to form-factors computations of correlation and response functions, where 
integrable techniques can provide, if not the exact result (which would involve resumming 
an infinite number of contributions) , controlled approximations of excellent accuracy, from 
the lowest energies through crossover regimes [^2| ||23|| . 

The physical origin of the accurate reproduction of scattering at non-zero energies 
relative to the Fermi surface at the Anderson model's symmetric point lies in a separation 
of scales. At this point there are two relevant scales in the problem: one is the Kondo 
temperature, Tjg, while the other is y/UT, a function of the Coulomb repulsion, U, and the 
resonance width of the dot level, V. At the symmetric point, <C y/UT. We exploit this 
fact to make our identification of integrable scattering states. But in turn this means that 
we expect errors in transport quantities involving scattering away from the Fermi surface 
of 0(T k /VUT). 

Bootstrapping from our success with the finite temperature linear response conduc- 
tance, we look at the non-equilibrium conductance near the symmetric point both in and 
out of a magnetic field at zero temperature. In order to compute this conductance we 
again employ a Landauer-Biittiker formalism akin to that employed in computing the out- 
of-equilibrium conductance of interacting quantum Hall edges p3}| . We imagine placing 
each lead at two differing chemical potentials, \i\ and \i<i- These differing voltages induce 
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different populations of free electrons in the leads. As with the finite temperature linear 
response problem, we recast these electrons in terms of our integrable scattering states. 
We then compute the (equilibrium) scattering amplitudes of scattering states in the leads. 
These scattering amplitudes then provide the probability for a state to tunnel from one 
lead to the other. Although the system is interacting, its integrability again implies the 
states scatter one-by-one. It is important to understand that this means of computa- 
tion introduces no additional error into the calculation. The sole source of uncertainty is 
found in whether the scattering states that we construct are entirely confined to a single 
lead. But because of the excellent agreement of the finite T linear response conductance 
with the previous NRG results, we expect this error to be similarly insignificant in our 
non-equilibrium computations. 

It is important to understand that this approach to the non-equilibrium physics differs 
from that used in |T2[ M [10] in a fundamental way. There the non-equilibrium conductance 



is expressed in terms of the non-equilibrium density of states of the impurity as determined 
from the correlator, lm(dd*). Here we have nothing direct to say about the non-equilibrium 
(or indeed, the equilibrium) behaviour of this quantity. 

We must stress this as the reader may be confused by the fact we do use the impurity 
momentum (which is in turn related to the Bethe ansatz impurity density of states as 
explained in Section 2) to compute the scattering amplitudes of excitations off the dot. 
This confusion may be heightened in that we employ the equilibrium Bethe ansatz density 
of states in computing the scattering matrices. It would thus seem legitimate to ask why 
we do not need to use a non-equilibrium impurity density of states in computing the out 
of equilibrium conductance [^5| . 

The answer lies in correctly understanding the basis of excitations by which we com- 
pute the conductance p6fl . We are able to use the equilibrium scattering phases as we 
employ the basis that is naturally present when the system is equilibrated. However be- 
cause of the integrability of the system, we can continue to employ this basis when we 
move the system out of equilibrium. These particles continue to scatter as they do in 
equilibrium. We note that if one were to compute out-of-equilibrium scattering matrices, 
one would find that they differ from their in-equilibrium counterparts by an overall phase 
alone. As transport quantities depend upon the absolute value of the scattering, this over- 
all phase would then have no effect |[26|| . While the application of finite voltage does not 
effect the scattering of the excitations, it does change their distribution in the leads. And 
indeed we must and do take this into account. 
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The rational behind this understanding has been tested beyond the various checks 
of these ideas applied by |24[] in their computations of conductances of quantum Hall 
edges. Generically the thermodynamics of an integrable field theory can be computed 
using thermodynamic Bethe ansatz which employs zero temperature S-matrices (and not 
finite temperature S-matrices as might again be naively expected - the finite temperature 
distribution might be thought to necessarily dress the scattering). With the thermody- 
namics, one can determine the finite temperature scaling behaviour of a field theory. If in 
the UV or IR limit, the theory flows to a conformal field theory, the finite scaling behaviour 
in these limits is independently determined by the central charge, c, of the theory. That 
the two computations always agree provides strong evidence we are handling the problem 
correctly. 

Turning to our non-equilibrium results, we find that they reproduce the expected gross 
features of the experimental differential conductance. When H = and we are near the 
symmetric point, the differential conductance is sharply peaked about its linear response 
value. The peak width is controlled by the scale, the Kondo temperature. We find 
that the peak is roughly symmetrical about V = in accordance with experiment [|l|. 

In a non-zero field, Meir and Wingreen |10[ predicted that the differential conductance 
would peak at eV = ±H. We find such peaks with our techniques although our peaks are 
found shifted to values of e\V\ notably smaller than \H\. Even in the limit of fields much 
larger than we do not find the peaks at \H\. This is again consistent with experiment 
0. We should not necessarily expect the peaks to occur exactly at eV = H as Meir 
and Wingreen's prediction is predicated in part upon a second order perturbative result. 
Moreover our construction of the scattering states suggests that in the particular case of 
the differential magneto-conductance, our results become exact in the limit of large applied 
fields. 



A portion of the results of this paper have been reported in |27fl • Here in this work we 
provide far greater detail on the nature of our computations. The paper is so organized as 
follows. In Section 2 we introduce the continuum version of the two lead Anderson model. 
The two lead Anderson model is integrable as is. However we first map it onto a one lead 
problem. If we were to explicitly solve the two lead problem, we would find nevertheless 
that we would be implicitly implementing the map to the one lead case. Having done this 
we review the Bethe ansatz for the one lead Anderson model together with the excitations 
necessary to form the ground state at zero temperature. The remaining portion of Section 
2 is devoted to identifying the excitations (both at and away from the Fermi surface) that 
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can be identified with scattering states and then computing their scattering amplitudes. 
We provide further details of the approximate nature of the scattering states so identified 
away from the Fermi surface. In the course of this discussion we demonstrate the Friedel 
sum rule. 

In Section 3 we explore the behaviour of the T = linear response conductance both 
in and out of a magnetic field. Because Wiegmann and Tsvelick E3] computed expressions 
for the occupancy of the dot, 7V|-/j., as a function of the gate voltage, Ed, and H in a 
variety of regimes, and the Friedel sum rule relates the electron scattering phase to this 
occupancy, we can derive closed form expressions for the linear response conductance in 
these same regimes. Outside these regimes we compute the occupancy numerically. Using 
these results, we show how the linear response conductance behaves as function of Ed- 

We plot the linear response conductance as a function of the gate voltage, Ed, in zero 
field in Figure 4. However the most interesting results of this section are found in our 
computations of the linear response conductance at finite H. In Figure 5 we plot the linear 
response conductance as function of Ed for a variety of values of H. As H is increased from 
zero we see that the linear response conductance is suppressed as a result of the destruction 
of the Kondo effect in a finite field. We also see the structure of the conductance peak 
evolve from its zero field value to that of free fermions. With increasing H, the full width 
at half maximum becomes narrower, decreasing from its zero field value of U to its free- 
fermion value of 2F. The peak height, in turn, decreases from its maximal value of 2e 2 /h 
to e 2 /h. And finally the location of the peak shifts from Ed = —U/2 to H/2 (we scale H 
here and throughout this paper such that gus = 1), appropriate to a spin up free-fermion 
with a field-induced, shifted chemical potential. These results are encoded in Figure 6. 

At the symmetric point U/2 = —Ed, we provide a simple closed form expression for 
the conductance (see Section 4.2). We then find that the conductance deviates from its 
unitary maximum for small fields via 

G = 2 e ^(l-^^ + (D(H*/T k )). (1.4) 

This deviation from the maximal conductance is quadratic, appropriate for the controlling 
H = strongly coupled Fermi liquid fixed point. Here Tj~, the Kondo temperature, is given 
in (3.16). 

In Section 4 we compute the finite temperature linear response conductance at the 
symmetric point, U/2 = —Ed- This requires recomputing the scattering of Section 2. At 
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finite temperature one must consider the thermal bath of all possible excitations. This 
highly non-trivial bath modifies the scattering. However doing so leads us to a highly 
pleasing result. We are able to reproduce Costi et al.'s NRG result for G as a function of T 
(see Figure 9). This is convincing evidence that we have correctly identified the scattering 
states (at least at the symmetric point of the dot). More specifically, we know the linear 
response conductance will again have a Fermi liquid form: 



Costi et al. demonstrated in perturbation theory that the constant, c, takes the value, 
c = 7r 4 /16 = 6.088. We compute in comparison c = 6.05 ± .1. Beyond the low temperature 
Fermi liquid regime, we emphasize we are able to describe accurately the conductance in 
the crossover regime T ~ T k . 

We also compare the scaling curve for the finite temperature linear response conduc- 
tance with the experimental data found in |l|. The comparison is plotted in Figure 10. We 
find the experimental measurements, even though taken away from the symmetric point 
of the dot (although still in its Kondo regime), agree well with our scaling curve. This 
suggests the Kondo regime of the dot exists over a wide range of gate voltages (i.e. Ed)- 

In Section 5 we move on to compute the non-equilibrium conductance at zero tem- 
perature. Again the scattering amplitudes need to be recomputed to take into account the 
change in the distributions of electrons in the leads induced by the finite bias. We also 
discuss subtleties with understanding how to think of a finite biased system in its one lead 
formulation. We then present results of the differential conductance both in and out of a 
magnetic field as discussed above. 

We are able to derive a number of simple closed form results for the out-of-equilibrium 
conductance. The differential conductance in zero field at the symmetric point is computed 
to be 



a remarkably simple result. We are also able to characterize the peak in the differential 
conductance that develops in the presence of a magnetic field. The quantity most discussed 
in the literature, the bias at which the peak occurs, is given by the expression 



G(T) = 2-(l-c^ + 0(TVT k )). 



e 2 , T 2 
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valid for H ^> and H <C yUT. Because of the logarithmic dependence of b upon H, the 
position of the peak approaches eV = \H\ extremely slowly. In addition to the location 
of the peak, we are also able to describe both the peak width and the height of the peak. 
The peak width is given by 

H ( 1 1 
eAV = — (tan -1 — — --tan" 1 — — -), (1. 



2ttV J -1 — 1/2 — b I~ 1 + l/2-b 

while the peak height is equal to 

gg (I-'-b) . 

°" h l 2 ^4(/-'-i.) + l J ' 1 ' 

We see that the height of the conductance peak approaches e 2 /h, one half the unitary limit, 
in the asymptotic limit of large H. Interestingly, the non-equilibrium Kondo model in the 



Toulouse limit [JUj also predicts a zero temperature differential magneto-conductance peak 
at eV = H which has a peak of e 2 /h, one-half the unitary limit. 



2. Basic Formalism 

2.1. Description of the System 

Pictured in Figure 2 is a sketch of the quantum dot connected to two leads. The 
Hamiltonian for this model in the continuum limit is given by 

/oo 
dx{-ic[ a (x)d x ci <T (x) + Vi5(x)[c\ a (x)d CJ + dlc la (x)}}+ 

« "°° (2.1) 

a 

where n a = d' a d a . Here ^ z is a sum over the two leads (/ = 1, 2). We have allowed for 
the possibility that the hopping matrix element, V;, differs between the leads as is typical 
in any experimental realization. Rather than treating the leads as half-lines with both left 
and right moving fermions, we represent the leads as 'unfolded' with fermions that are 
solely right-moving. Fermions in either lead that are incident upon the dot are considered 
to lie in the region, x < 0, while those traveling away from the dot in either lead are found 
with x > 0. We represent this in Figure 2 by drawing the leads as elongated arcs. 
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Figure 2: A sketch of two leads attached to a quantum dot. 



We stress that there are no interactions in the leads (as opposed to the calculations in 



24] for instance): the non-trivial physics of the problem arises solely from the strongly 
interacting dynamics of the dot. 

It will be advantageous to reformulate this problem as a one-lead Anderson model 
(i.e. a single lead model). To do so, we introduce even/odd electrons: 

1 



c e/o = -^^(V 1/2C1 ± V 2/1 c 2 ). (2.2) 
Recasting 7i in this new basis, the odd electron, c c , decouples and we are left with 



H = J2 I dx{-icUx)d x c ea (x) + (V? + V 2 2 ) 1 / 2 5(x) [cl <7 (x)d a + dtc ea (x)]}+ 

(2.3) 

a 

We have thus reduced the problem to that solved using Bethe ansatz in a series of papers 
by Kawakami and Okiji |0| and Filyov, Wiegmann and Tsvelick [|16 . 

With this reformulation of the model, we have to address the question of computing 
scattering amplitudes of electronic excitations in the original two lead problem. Naively it 
would seem we can do this. Let T(e) and R(s) be defined by: 

T(e) : transmission amplitude for an electronic excitation of energy, e, 
to scatter from lead 1 to lead 2 (or 2 to 1), 

(2.4) 

-R(e) : reflection amplitude for an electronic excitation of energy, e, 
to scatter from lead 1 to lead 1 (or 2 to 2). 

Assuming the amplitudes behave linearly under the even/odd map ( |2.2|) , we have (for 
V 1 = V 2 ) 

jSAe) =R( £ ) + T(e); 

(2.5) 

e iS (e) = i = i2( e )_r(e). 
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We then conclude that 



e iS e (e) _ j 

T(e) = ; 

a A ( 2 - 6 ) 

= / - , 

govern scattering in the two lead problem. However not all electronic excitations behave 
linearly under the map (|2.2| ) and so things are not always this simple. We will consider 
this in more detail in Section 2.4 and in particular Section 2.5. 

As noted, in writing the above two equations, we have assumed V\ = V%. If V\ ^ V2, 
the transmission amplitude is scaled by the factor, 

2V X V 2 



2 • 



As this is a constant factor, it always possible to rescale results to take into account an 
asymmetry in the dot-lead couplings. As such, we will assume throughout the paper that 
V x = V 2 . 

2.2. Bethe Ansatz Solution of the One-Lead Anderson Model 

In order to understand the scattering between the two leads we will rely upon aspects 
of the one- lead Bethe ansatz solution found in [ 15 1 and |16|] . As such we summarize briefly 



the results of this work. Applying the Bethe ansatz yields a set of quantization conditions 
describing a finite number of bare excitations in the system: 



e 



n^.)- Aa -z/2' 



where 



1^- A Q - g(kj) + i/2 = _ tt A Q - A/3 + i 
f = \ Aq - g(k 3 ) - i/2 ~ 11 A a - A^ - i ■ 



6(k) = -2 tan" 1 
9{k) 



(2.7) 



(k - e d ) 

(k-e d -U/2)\ (2.8) 
2UT 

r = (v 1 2 + vi). 

As in all problems with an SU(2) symmetry, there are two types of excitations: charge 
(with rapidities, k) and spin (with rapidities, A). Here iV is the total number of particles 
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in the system, and M marks out the spin projection of the system, 2S Z = N — 2M (in zero 
magnetic field, M = N/2). 

When Ed > —U/2, the ground state of the system consists of the following set of 
excitations: 



i) N - 2M real kjs; 

ii) M real A^s; 

iii) associated with each of the M A^s are 

two complex k's, k±, described by (2-9) 
g{k%) = g(x(X a ) q= iy(X a )) = A a ± i/2; 
x{\) = U/2 +e d - VUf(X + (A 2 + l/4) 1 /2)i/2. 
y(X) = VUf(-X + (A 2 + 1/4) 1 / 2 ) 1 / 2 . 

Although only valid for Ed > —U/2, we can also understand the case, Ed < —U/2, through 
a particle hole transformation. If we take the continuum limit of ( |2.9|) , we no longer deal 
with discrete values of A and k, but rather go over to smooth distributions, p(k) for the 
real fej's, and cr(A) for the A Q 's and their associated fe^'s. To derive these distributions, 
we first take the log of (|2.7[): 



M 



kjL + 8(kj) = 2kN 3 - 0i(9(kj) - Xp); 



[2.10) 



13=1 

M N-2M 

2nJ a + ^e 2 (x a -Xp)+ ifr« - gfc)) 

0=1 3=1 

= -2Lx(A a ) - 2Re5(x(A Q ) + iy(X a )); 
2 

eJx) = 2tan _1 (-x) + w. 

n 

Nj and J a are the quantum numbers of the charge and spin excitations respectively. 
Taking the thermodynamic limit (i.e. N, M, L — > oo with N/ L and M/ L finite) , followed 
by derivatives of the above, gives 



1 A(k) f Q 
P(k) = — + -jr 1 + g'{k) / dX a x {g{k) - A)cr(A); 
2-k L J Q 

a(X) = + ^ - [ Q dX'a 2 (X' - X)a(X') - 

* L Jq 



B 



(2.11) 



D 



dk ai(A — g(k))p(k), 
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where 

A(k) = ^-d k 5(k); 
LlX 

A(A) = ——d\Red(x(X) + iy(X)); (2.12) 

7T 

a n{z) 7; d x 9 n \x) -. ~ ■ - lyT- 

An 7T [n z + 4x^) 

Various limits appear in the above equations for the distributions. — D marks the lower 
allowed limit of the /c's while Q marks out the bandwidth of the A's. As each A has a pair 
of complex /c's, its associated energy is 2x(A). We thus determine Q by 

x(Q) = -D. (2.13) 

Often it will be possible to replace Q and — D by 00 and —00. B and Q on the other hand 
give the spin and charge Fermi surfaces. They are determined by the constraints 



N -2M 



B 

dk p(k); 

-d 

o (2 - 14) 

m r Q 

— = / dX cr(A). 

2.3. Determination of the Scattering Phase at the Fermi Surface: The Friedel Sum Rule 

In this section we examine the relationship between the scattering phase of electrons 
S e (e) at the Fermi surface and the number of electrons on the dot and so verify the Friedel 
sum rule. 

To determine S e (e), we employ an energetics argument of the sort used by N. Andrei 



in the computation of the magnetoresistance in the Kondo model |T?]. Imagine adding an 
electron to the system. Through periodic boundary conditions, its momentum is quan- 
tized, p = 2-Kn/L. If the dot was absent, the quantization condition would be determined 
solely by the conditions in the bulk of the system and we would write, ptmik = Aim/L. 
Upon including the dot, this bulk momentum is shifted by a term scaling as 1/L. The 
quantization condition is then rewritten as 

2tt^ S e (e) 
P=-^=Pbuik + ^, (2.15) 

where L is the system's length. The coefficient of the \ term is identified with the scattering 
phase of the electron off the dot. 
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As we are interested in expressing 5 e in terms of the number of electrons on the dot, 

it is useful to separate out from p(k) and cr(A) the impurity contribution to the density of 

states. We thus write ^ 

p{k) = pbuik(k) + yPimp(k); 

\ (2-16) 
cr(A) = abuik(A) + — cr imp (A). 

Pbuik/o"buik represent the bulk contribution to the densities while pimp/cimp determine the 
number of electrons of definite spin, n^/ridi, sitting on the dot. From ( |2.14j ) we have 

i dT = dX cr imp (A) + / dk Pim P (k); 

\ J ~ D (2.17) 

Q 



r^l 

n di = dX cr imp (A). 

JQ 



These relations will be key in verifying the Friedel sum rule. Substituting ( |2.16| ) into ( |2.11|) , 
we obtain separate equations for pbuik/pim P and <7buik/cimp: 

1 

— +g'(k) dX ax(g(k) - A)cr bu i k (A); 
27r Jo 



Pbulk(^) 



Cbuik(A) 
and 



x'(X) 



n 



rQ rB 

/ dX'a 2 (X' - A)cTbuik(A') - / dk ai(A - g(k))p hu i k (k) 
Jq J-d 



(2.18) 



Pimp(k) = A(k) + g'(k) / dX ai(g(k) - A)a imp (A); 

Jq 

rQ r B 
0i m p(A) = A(A) - / dX'a 2 (X' - A)cr imp (A') - / dk ai(A - g(k))p iuip (k). 
Jq J-d 



(2.19) 



In Appendix A we give alternative forms to the above equations governing the density 
functionals. These alternatives are far more amenable to numerical analysis and in practice, 
the ones used in solving for the densities. 

Having obtained the equations governing the impurity densities of state, we now focus 
on the scattering phase itself. From (|2.10|) we can read off the bulk momentum of a 
charge/spin excitation with quantum number N/ J to be 

p(k) = = k+ dX a hulk {X)e t {g(k) - A); 

Jq 



L 



Q 

B 



p(X) = = 2x(X) + [ dX'a hulk (X')e 2 (X -X')+ [ dk p hu i k (k)d x (X - g(k)). 

L Jq J-d 

(2.20) 
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We assume here tan -1 in varies from n/2 to n/2 so ensuring a simple relationship 
between the momentum and the energy functionals to be derived in the next subsection. 

The impurity contribution to the momentum for each type of excitation can be simi- 
larly determined to be 

pQ 

Pim P (k) = 6(k) + / dX a imp (X)(6 1 (g(k) - A) - 2tt); 

JQ 



Q 

•Q 



Pimp (X) = 2Re5(x(X) + iy(\)) + [ dX'a irnp (X')(6 2 (X - A') - 2tt) ( 2 - 21 ) 

JQ 

+ [ dk p imp (k)(6 1 (X-g(k))-2n). 

J-D 

Here we have chosen a different range for tan -1 in for describing the impurity mo- 
mentum. Shifting back to the original range leads then to the appearance of the 27r's. 
This choice is governed by our ultimate desire to give the scattering phases in terms of the 
impurity momentum. In particular we want p\ mp {k — > — oo) = Pi mp (X — > oo) = 0. 

According to ( |2.15| ), we identify Pi mp (&) with the scattering phase of a charge excita- 
tion and pimp (A) with the scattering phase of a spin excitation. By differentiating these 
expressions and comparing to (|2.19|) , we obtain the relations 

dkPimp(k) = 2np imp (k); 

(2.22) 

d\Pi m p(X) = -27ra imp (A). 

Again we have relations crucial to verifying the Friedel sum rule. 

In order to determine the scattering phase of an electron (as opposed to a spin or charge 
excitation), we must must specify how to glue together a spin and a charge excitation to 
form the electron. The situation is analogous to adding a single particle excitation in the 
attractive Hubbard model Adding a single spin f electron to the system demands 

that we add a real k > B (charge) excitation. In doing so we create a hole at A > Q in the 
spin distribution as the number of the available slots in the spin distribution is determined 
by the number of electrons in the system. Adding an electron to the system thus opens 
up an additional slot in the A-distribution. 

The electron scattering phase off the impurity is then the difference of the right- 
moving k-impurity momentum, Pi mp (/c), and the left-moving A-hole impurity momentum 
-Pimp (A): 

S l = Pimp = Pimp(k) + Pimp (A); 

rh f Q (2.23) 

= 2tt / dk'p imp (k') + 2n / d\'(T imp (\'), 

J-D J\ 
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where we have used ( j2.22|) in writing the last line. If the excitations are added/removed at 
the Fermi surfaces, i.e. k = B, A = Q, we obtain the Friedel sum rule for spin up electrons, 



8\ = 2tt / dkp imp (k) + 2tv / dXa iuip (X) = 2im dh (2.24) 

J-D JQ 



where (|2.17|) has been used. The total energy of this excitation is e(k = B) + e(X = Q) = 0, 
as it should be. 

To determine the scattering of a spin down electron we employ particle-hole symmetry. 
A particle-hole transformation is implemented via 

-» ci(-k); 
c\(k) -> c T (-fc); 

d\^d i; (2.25) 

d\ -> d T ; 

£d -»• -U - e d . 

Consequently the scattering phase of a spin { hole is related to that of a spin ] electron 
via 

Sl(-U-e d ) = Sl(e d ). (2.26) 
The phase of this excitation is then 

S l ho (-U - e d ) = 2tt [ Q dX'a imp (X') + 2n [ dk'p imp (k') = 2nn d] ie d ). (2.27) 

JX J-D 

where the last equality holds if we take the hole to be at the Fermi surface, A = Q and 
k = B. As n d] {e d ) = 1 - n dl (-U - e d ), we have 

S l ho (-U - e d ) mod 2tt = -n dl {-U - e d ). (2.28) 

At the Fermi surface, hole and electron scattering are identical (up to a sign) and so we 
verify the Friedel sum rule for spin down electrons. 

The reader may be puzzled why we rely on a particle-hole transformation in computing 
the scattering amplitude for spin-down electrons. Although it would be desirable to do 
this computation directly, it does not seem to be possible. To construct a spin { electron 
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at the Fermi surface, it is natural to remove a k = B excitation while adding a A = Q 
excitation. The corresponding scattering phase is then given by 



5 l e = pi 



mp 
Pimp 

(k) + 

Pimp (A) 

f B ,Q (2.29) 

2tt / dkp irnp (k) + 2n d\'a imp (X') 

J-D JQ 



= 27m dT . 

But this is obviously not what we want - a manifest violation of the Friedel sum rule. 
Rather by comparing ( |2.29| ) with ( |2.27| ), the scattering indicates we have constructed a 
spin I electron not at Ed but at the particle-hole conjugate point, —U — Ed- Why this is so 
it is not entirely clear. However one can notice that the k-excitations are not only charge 
excitations, but are in some sense unbound spin j electrons (the number of k-excitations 
is directly proportional to the magnetization of the system). So in removing a k-excitation 
to form the spin j electron, we are in some sense creating a spin j hole. And a spin f hole 
at chemical potential, e^, will scatter as a spin j electron at —U — Ed- 

This entire discussion has concerned itself with proving the Friedel sum rule for the 
one-lead Anderson model. However we can argue that it applies, appropriately revised, to 
the two-lead model. More precisely, we can argue that the excitations at the Fermi surface 
behave linearly under the map ( |2.2| ) and so have two-lead scattering amplitudes given by 
( |2.6|) . This will be detailed in the following two sections. 

In appendix B we give an alternate derivation of the scattering phase that focuses 
upon the impurity energy of an excitation as opposed to its impurity momentum. In doing 
so, we elucidate subtleties not explicitly discussed in f28[j . We also give a third derivation 



of the scattering phase in Appendix C by directly considering the dressing of the bare 
scattering. 

2-4- Excitations Away From the Fermi Surface in the Kondo Regime 

In the previous section we were mainly concerned with scattering at the Fermi surface. 
However as made clear by taking k ^ B, A 7^ Q, we can look at scattering above the Fermi 
surface. 

It is tempting to ask first whether the non interacting electrons in the lead can still 
be described in this formalism (by electrons we mean the standard plane wave excitations 
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of appropriate spin and charge). Here it is useful to recall some well known results from 
many-body theory: Langreth, in verifying the Friedel sum rule for H = ||3Q|| . computed 
the ratio of the elastic inverse lifetime, tT 1 , of a plane wave mode to that of its total 
inverse lifetime, r _1 , finding 

r-He) ~2bnE(eY 1 J 

where is the self-energy for the dot electron Green's function. At the Fermi surface, 

lm£(e = 0) = ~ (2.31) 

and there are no inelastic processes. However away from the Fermi surface 

p 

ImE(e) = - + ce 2 , c> 0, (2.32) 

and rj 1 < r _1 , so electrons with energies above the Fermi surface do not scatter elastically. 

On the other hand, the simple excitations we construct within the integrable descrip- 
tion by gluing spin and charge excitations will necessarily scatter elastically: beyond the 
Fermi surface, they cannot be the free electrons one would initially like to describe. In and 
of itself, this does not matter as all we are interested in at the end is charge transport, ir- 
respective of what kind of objects actually do carry this charge. A similar situation occurs 
in the fractional quantum Hall effect |24[], where the integrability approach uses quasi- 
particles which are neither electrons nor Laughlin quasi-particles. This approach merely 
provides a more convenient basis for the space of excitations, chosen such that scattering 
at the impurity is as simple as possible. However we do have to concern ourselves with 
rewriting the original free electrons in terms of the integrable scattering basis. As indicated 
in the introduction we are not able to provide an answer to this problem in its entirety. 

Our approach will then be to build excitations which are "electronic" , that is carry the 
same quantum numbers as electrons, but scatter simply (i.e. elastically) at the impurity - 
they will also scatter in a simple, factorized way among themselves, although their S-matrix 
is non trivial (it is not S = — 1 anymore). One can certainly think of these excitations as 
dressed electrons. 

This being understood, another difficulty remains: the potential parameter space of 
electronic excitations, i.e. (k,X), is two dimensional (provided we neglect other solutions 
of the Bethe ansatz, see Section 2.6), whereas we naturally want the space to be one 
dimensional. For the moment, we can only make the necessary dimensional reduction 
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when we are in the Kondo regime of the Anderson model. The first step in doing so is to 
determine the energy-momentum of an excitation labelled by (k,X). 

We already know the momentum of the excitations from (|2.20| ). We thus must only 
compute the energies. To facilitate the calculation of excitation energies, it is useful to 
decompose p(k) and <r(A) into particle and hole densities: 



energy. We can write this variation in two ways: one in terms of the bare energies and one 
in terms of new functions, e (k) and £ ± (A), governing the dressed energies: 



5E = L dk {e+(k)5p p (k) - e-(k)5p h (k)} + L dX {e + (X)5a p (X) - e~ (X)5a h (X)} 



p p/h (k)=e(±B T k)p(k); 
a p/h (X)=6( T Q±X)a(X). 



(2.33) 



Now imagine varying p p /h and cr p /h and asking what is the corresponding variation in the 




(2.34) 



The variations on 8p p /h and Sa p /h are not independent. From ( |2.11|) we see 




(2.35) 



Substituting ( |2.35| ) into ( |2.34j ), we obtain 



e+(k) +e~{k) = k- — - J dX e~(A)ai(A - g(k)); 
e+(X)+e-(X) = 2x(X)- [ dX'e~ (X')a 2 (X' - A) 



(2.36) 




e ± (A) and e ± (k) are characterized by 



e+(X) = 9(Q - X)(e+(X) + e"(A)) > 0; 
£ -(A) = 9(X - Q)(e+(X) + £"(A)) < 0; 
e+(fc) = 0{k - B)(e + (k) + e~(k)) > 0; 
e~(k) = 6{B - k)(e+(k) + £~(k)) < 0. 



(2.37) 
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The functions, e = e + + e~ , are continuous and monotonic. e have been defined such 
that e + (k/X) is the cost of adding an excitation at k/X while —e~(k/X) is the energy 
needed to create a hole at k/X. Again in Appendix A, we give alternative forms to the 
above equations governing the energy functionals which are more amenable to numerical 
analysis. 

Having determined the energy of the excitation, we can easily relate it to its corre- 
sponding momentum. We consider the case of H = first. Comparing (|2.20|) and ( |2.36|) , 
and using dke(k) = 27rpbuik(&) and d\e(X) = — 27ro"b u ik(A), we see 

Pbuik(k) = e(k); 

(2.38) 

Pbuik(A) = e(A), 
where pbulk is the portion of momentum not scaling as 1/L. 



k 




Figure 3: A cartoon of the parameter space describing electronic excitations. The 
drawing supposes that Sq = < S\ < e<i < £3 < £4. Each curve represents a set of 
the excitations that share the same energy e. Only in the case, Eq = 0, i.e. when 
we are at the Fermi surface, is the pair (k, X) uniquely specified. The dashed line 
marks out the ansatz we employ in the Kondo regime. 
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With this in hand, we can parameterize the scattering phases of electronic excitations 
away from the Fermi surface. Suppose we want to characterize a spin f electron with 
energy, e e [. (This is sufficiently general for H = as we know spin j. electrons will scatter 
identically.) The possible k and A forming this excitation must satisfy 

e(k) - e(A) = e ei . (2.39) 



Given (|2.38| ), this choice automatically satisfies e e i = p e i (up to 1/L corrections). 

This parameterization leaves an unresolved issue. It does not in general uniquely 
specify a particular pair (/c, A), crucial if we are to actually compute quantities involving 
information away from the Fermi surface. We have schematically illustrated the degen- 
eracy of choices in Figure 3: as the energy is increased the multiplicity of pairs (k, A) 
correspondingly increases. 

In certain cases, however, the specification is unique. At the Fermi surface, the degen- 
eracy of pairs is lifted. This already illustrated in Figure 3. However there is another case 
not captured by the cartoon in Figure 3 where the specification is unique. We have only 
a single possible pair (k, A) for each energy in the case of spin j hole scattering in finite 
magnetic field at the symmetric point of the model. The reason behind the reduction of 
the parameter space for this case will be made clear in what follows. 

For the other cases where the choice is not unique, the question becomes on what 
operative principle do we reduce the space. The key that we have identified to reducing 
the parameter space is determining how the excitations behave under the map fl2.2j) . Only 
if they behave linearly under the map are they of use for it is only then that we can 
compute their scattering amplitudes in the two- lead picture via ( PHI) . 

Although we have refined the question, we cannot in general determine whether a 
given excitation unfolds linearly in the case when there are multiple possible pairs (k, A). 
We can however make some progress when we are in the Kondo regime of the Anderson 
model (i.e. U + 2e^ ~ 0). In this regime we expect the scattering phase to vary on the 
scale of the Kondo temperature, T^. The electron scattering phase is determined by pi mp 
and fTi mp , the two impurity densities. Of the two, only pj mp varies on scales on the order 
of Tfc. In contrast, a; mp is controlled by the much larger scale yUT. Thus in computing 
electronic scattering phases away from the Fermi surface at zero temperature, it is natural 
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to keep A = Q, its Fermi surface value, and vary k. Specifically, to describe an electron of 
energy, e e i, we chose (k, A) such that 

k particle, e(k) = e e i; 

(2.40) 

A hole at A = Q. 

With this ansatz, we then have restricted the two dimensional phase space, (A, k), of 
potential excitations carrying the quantum numbers of an electron to an one dimensional 
subspace. Hence the scattering phases of electrons of energy, e e i, above the Fermi surface 
at H = are given by 

5l(e el ) = 5i(e el ) = 2n ! dk'p imp (k') + 2n C dX a imp (A), e(k) = e el . (2.41) 

J-D JQ 

When if ^ 0, we still have a simple relation between the energy and momentum, i.e., 
we have 

«(*)-!<*)- fi (2 42) 

e(\) = p(A). 

Hence with H ^ we are still faced with an oversized parameter space. But we conjecture 
similar relations to those in (|2.40|) hold in constructing the electronic spin f excitations: 



| electron : k particle, s(k) = e e f, 

(2.43) 

A hole at A = Q. 



The scattering phase of this excitation is accordingly 

6l(e el ) =2n [ dk'p imp (k') + 2n f Q dX a imp (A), e(k) = e el . (2.44) 

J-D JQ 

With H 7^ and consequently e~(k) not identically zero, we can construct spin f hole 
states by removing a /c-state and a A-hole. The scattering phase of spin j holes is then 
equal to 

S ho( £ ho > 0) = Pimp(k) +p imp (X); 

rk rQ (2.45) 

= 2n dk'p iiap (k') + 2n / dX'a iiap (X'), e(k) = -e ho . 

J-D JQ 

For this particular excitation we do not need the scattering ansatz. If we are to remove 
a A-hole we must do it for A < Q. However at the symmetric point Q = — oo and so 
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the choice is unique. This fact will allow us to conclude that in the large field limit, our 
computation of the differential magneto-conductance becomes exact. We also point out 
that as e{k) is bounded below, i.e. e(k = —D) = —H, we are limited in the energy range 
([0, H\) in which we can construct spin f holes. 

With H 7^ 0, we must compute the scattering of spin { objects separately. To do so 
we again employ a particle-hole transformation. We so obtain 

S^(e e i, -U - Ed) = Si {E ho = £ eZ , e d ); 

i T (2 - 46) 

5ho( £ ho, —U - Ed) = 5l(e e l = Eho, Ed). 

Unlike scattering at the Fermi surface where we could infer behaviour at Ed > —U/2 from 
behaviour at Ed < —U/2, we cannot do so for scattering away from the Fermi surface. Thus 
when we need to compute explicitly the scattering of spin j objects, say in computing the 
magneto-conductance out of equilibrium, we will be restricted to the symmetric point 
E d = -U/2. 

2.5. Returning to the Two-Lead Problem 

In this section we explore in more depth the map between the one and two lead 
models and its attendant problems. To review the map in more formal terms call E Bj0 the 
integrable excitations in the even and odd leads. We then describe the factorized scattering 
by the relations 



E e D = e lSE DE e 
E Q D = e l5 °DE , 



(2.47) 



where D is a formal symbol representing the impurity. Again the phase, 5 e , is non trivial, 
while S = 0. Under the map ( j2.2p , integrable excitations in the two lead picture are given 
by 

£1,2 = E e ± E . (2.48) 
Scattering of an excitation in the first lead is then described by 

E l D = RDE 1 +TDE 2 ; 

or 

(E e + E )D = RD(E e + E ) + TD(E e - E Q ), (2.49) 

where consistency with ( |2.47|) demands that the transmission and reflection amplitudes, R 
and T, satisfy R+T = e lSe and R — T = e lS ° . An implicit assumption in this determination 
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of the scattering in the two lead picture is that the superposition E Q + E e of an electronic 
excitation in the even sector and an electronic excitation in the odd sector carries unit 
charge in lead 1 and no charge in lead 2. For an arbitrary fermionic excitation in the even 
and odd leads this will not be the case. For example, imagine a electronic excitation in 
the even lead, that if decomposed into a plane wave basis of free electrons, consists in part 
of particle-hole excitations: 

E e = akcjk jFermi sea) + ^ a kk p k h c ek c e k p c ek h |Fermi sea) H , (2.50) 

k k^kp^kji 

where here c e fe is a plane wave electron in the even lead with wave vector k. The linear 
combination E e + E a , where the excitation, E a , is arrived at from ( [2.50|) through c e — >• c G , 
does then not carry unit charge in lead 1. Rather it carries indefinite charge in both leads. 
In this case the excitation does not transform between the two pictures as indicated by 
( p.48|) and its scattering cannot be expected to be described by (|2.49|) . If E e was strictly 
a linear combination of terms of the form c^lFermi sea), this problem would not surface. 

Thus in order to exploit the map between the two pictures, we must limit our self to 
excitations that behave linearly under the map as described in (|2.48 ). There are then 



two questions to be answered. Do such excitations in general exist? And if they do exist, 
are they sufficient for our purposes, the description of transport properties. We have two 
arguments that excitations with scattering described by (|2.48| ) do exist. More precisely we 
have two arguments giving that excitations falling along some line in the two dimensional 
(k, A) parameter space have such scattering. Moreover, we argue that the scattering of 
such excitations is sufficient to determine transport properties. 

The first argument relies upon the transformation properties of the Fermi field, c(x). 
Recall that in order to implement the map between the one and two lead pictures, it is 
c(x) that is transformed, i.e. 



c e (x) -> -y=(ci(x) ± c 2 (x)). 
Thus any integrable excitation E e that has a finite matrix element with c e , i.e. 

( Ce \E e ) ^0, 

must also behave linearly under the map from one to two leads. To see this more explicitly 
imagine making a mode expansion of the field, c e (x), in terms of the integrable excitations. 
As E e couples to c e , E e must appear in this expansion, 

c e (x) = J2 a re iPrX E er + ---, (2.51) 
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where E e is one of the E er J s. As c e and E e are linearly related, they must share the same 
transformation properties. This guarantees that any excitation appearing in the above 
mode expansion will have scattering described by ([2.49 ). 

But we can say more on the basis of the properties of c e (x). Because the underlying 
model is essentially free, we know the single particle spectral function of the model will be 
given by 



Thus for any given energy E = p, we know that some integrable excitation with this 
energy and momentum must appear in the mode expansion (|2.50| ). In terms of the two 
dimensional parameter space, (k, A), this implies that there is at least one line in this space 
describing excitations transforming as (|2.48| ) and scattering as ( |2.49|) . 

The second argument for the existence of this line in the (k, A) parameter space relies 
upon combining the properties of the low energy sector of the theory with the equivalence of 
the integrable excitation we have constructed at the Fermi surface with the corresponding 
plane wave electron excitation. Given our reliance on this equivalence, it deserves further 
exploration. 

This equivalence is, of course, strongly suggested by our ability to reproduce the 
Friedel sum rule and the fact Langreth [|30|] demonstrated that plane wave electrons at 
the Fermi surface scatter elastically, the hallmark of integrable excitations. Neverthe- 
less, the statement that the integrable excitation coincides with a plane wave electron 
needs further clarification. If we denote the wave function of the integrable excitation as 
4>int(x, x\, . . . , xn), where x is the coordinate of the excitation and the Xi are the coor- 
dinates of the electrons in the Fermi sea, and ^free e \.{x, ■ ■ ■ , %n) as the many-body 
wave function of the corresponding plane wave electron plus Fermi sea, we know that the 
orthogonality catastrophe implies the matrix element, (int|free el.), equals 



and so vanishes in the thermodynamic limit. Thus it would seem that that in fact that 
the two excitations do not coincide. 

However we are not interested in matrix elements involving full eigenstates of the 
Hamiltonian, but matrix elements involving asymptotic scattering states defined far from 
the impurity. These are the states of concern in applying a Landauer-Buttiker formalism. 
With such states, we would evaluate the above matrix elements by restricting x, Xi <C or 



(c e c e (0))(E,p)cx8(E-p). 




oo 



dxdxi ■ ■ -dx N il;* nt (x,xi, . . . , x N )ip{ ree ei.(^,^i, • • -,xn) = ^(1/L), 



— oo 
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x, Xi 3> 0, depending on whether the state is in-going or out-going. With such a restriction, 
the orthogonality catastrophe does not apply and 



(int|free el.) = 1 + 0(1/L). 

In this sense the excitations coincide. 

With this equivalence so understood, the two excitations, the integrable excitation 
and the plane wave share the same transformation property (|2.48 ) under the map. We 



now exploit this fact by combining it with the behaviour of the low energy sector of the 
theory. In this sector we can take a scaling limit and obtain a relativistic theory invariant 
under Lorentz transformations. Under such transformations, we can imagine boosting the 
integrable excitations at the Fermi surface, obtaining in the process an excitation with 
finite energy and momentum. However the transformation properties of the excitation 
cannot be altered by the boost. As such the boosted excitation will still transform via 
( ET48|) under the map. This again implies that there is a line in (fc, A) parameter space 
describing excitations transforming in the desired fashion. In this case moving along this 
line amounts to making a Lorentz boost. 

Having argued that there do exist excitations transforming as (|2.48|) , we now have 
to address whether the existence of such excitations is sufficient for our computations of 
transport properties. We can answer in the affirmative. To compute any given transport 
quantity in the Landauer-Buttiker approach, we need to sum up transmission amplitudes 
over some given energy range. For example, if we were to compute the zero temperature 
out-of-equilibrium conductance, this energy range would be determined by the difference 
of chemical potentials in the two leads. Now imagine looking at a particular infinitesimal 
energy interval within this range. As we know the density of states of the free electron in 
the lead, and we know that the interactions in the problem do not affect this particular 
quantity, we know precisely how much charge lies in this interval. Now we are able to 
construct an integrable state that transforms via (|2.48|) with an energy in this interval 
and with the same density of states as the free particles. Thus our integrable state com- 
pletely exhausts the charge lying within this infinitesimal interval. Given that we are able 
to compute its transmission amplitude ( J2.49[ ), we can compute the contribution of this 



infinitesimal energy interval to the transport quantity. 

As a corollary to this, the manifold other integrable states arising from the Bethe 
ansatz equations are then not needed for the computation of transport properties. We 
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do not need to account for the (k, A) states that do not transform as (|2.48| ). We also do 
not need to worry about states consisting of (k, A) excitations together with particle-hole 
excitations of k and A, or indeed excitations involving more complicated string solutions 
of the Bethe ansatz equations. The inclusion of such states in the computation of any 
transport quantity would amount to a double counting, given that the line of (k, A) states 
transforming as ( |2.48| ) completely exhausts the density of states of free electrons in the 
leads. 

Given all of this, we still must stress that our computation of scattering amplitudes 
away from the Fermi surface is in general only approximate. Although we understand that 
there exists a line of excitations in the (k, A) parameter space for which we understand 
and can compute scattering, we do not know which line. Rather, at the symmetric point 
of the model we have only an ansatz of how this line cuts through parameter space. 
However we again stress that this ansatz is supported by the nature of the two scales in 
the Kondo regime, 2*. and y/UT. Moreover our ansatz appears to be extremely good given 
its agreement with Costi et al.'s NRG results. 

Fortuitously there is one case where this ansatz is exact: the description of spin f 
holes. There, the (k, A) parameter space is one-dimensional from the start and no ansatz 
is needed. We point out that the scattering of such holes provides by far and away the 
main contribution to the differential magneto-conductance at large fields, H. As such we 
expect the differential magneto-conductance to be exact at asymptotically large fields. 

In order to determine for all cases how the line of linearly transforming excitations cuts 
through the (k, A) parameter space we would have to have complete control of the change 
of basis between the free electrons and the integrable excitations. There must presumably 
exist a complete set of such excitations providing a proper change of basis of the form 

\E) = J2cn,F 1 ...F k \F 1 ,...,F k ). (2.52) 

Here, the notation is highly symbolic: E stands for any possible integrable excitation, while 
Fi stands for free electrons with some spin and energy with the sum over the number, k, the 
types, and the energies of such excitations. Presumably, there are some complex selection 
rules making it possible to match the number of degrees of freedom on the left and on 
the right of ( |2.52|) . At this point however, an exact understanding of ( 2.52 ) is beyond our 
reach. 
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3. Linear Response Conductance at T = 

In this section we discuss the linear response conductance both in and out of a magnetic 
field at zero temperature. We of course note that none of Section 2 is necessary to compute 
the linear response conductance at T = 0. Although we have demonstrated the Friedel 
sum rule using integrability, all we need for this quantity is the occupancy of the dot as a 
function of various parameters, something available from the original Bethe ansatz work on 
the one-lead Anderson model. However, the behaviour of the linear response conductance 
as predicted by the Bethe ansatz has never been adequately explored, particularly in 
the case with a magnetic field. Indeed in the original work of Ng and Lee showing 
that the Friedel sum rule could be applied to quantum dots, they employ a Hartree-Fock 
approximation in estimating the dot occupancy and so obtain some qualitatively incorrect 
predictions as to the behaviour of the conductance. 

Generally, the linear response conductance equals 

G=^(|T T | 2 (e = 0) + |Tj 2 (e = 0)), (3.1) 

where Ty^(e = 0) is the scattering amplitude at the Fermi surface: 

\Tyi(e = 0)| 2 = sin 2 (nn dyi ). (3.2) 

In the case with H = 0, the number of electrons on the dot as a function of Ed, the gate 
voltage, can be computed exactly, as has been done by [p8[ . When H ^ 0, the equations 



become more difficult to analyze and in general only numerical solutions are available. 
However at the symmetric point, it is again possible to compute in closed form the number 
of spin t / I electrons on the dot [p8| , and so arrive at an analytic expression for G. We 
first consider the case with H = 0. 

3.1. H = Linear Response Conductance 
In this case the Friedel sum rule tells us 

\rp ,2 ■ 2 = K / Q q\ 

l T T/J = sm (^^ )> I 3 - 3 ) 

where the phase, <5f/j_, is equal to 



8yi = 2im d yi. 
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(3.4) 



The number of electrons, n^/i, on the dot when H = simplifies to 



n dyi = / cL\cr imp (A). (3.5) 
Jq 

Cimp(A) in turn is given by ( |2.19|) with the charge Fermi surface, £?, set to the bottom of 
the band, — D, 

a imp (A) = A(A) - [ Q d\'a 2 (X' - A)a imp (A'). (3.6) 
Jq 

The Fermi surface, Q, of the spin excitations is determined by the equations 



N 

T 



Q 

dX a(A); 

Q 



Q (3 - 7) 
a (\) = -±¥±_ f dX'a 2 (X' - X)a(X'). 

* Jq 

These equations are solved explicitly over most of the relevant parameter range in |28| using 
a Wiener-Hopf technique. 

The solution breaks down into three cases according to the value of Q describing the 
Fermi surface: 

case i: If we are close to the symmetric point (U/2 + Ed ^ y/(UT)) then Q <C (at the 
symmetric point, Q = — oo) and we have 

i i 00 ( —A\ n r°° 
n dm = \-^=Y. W^S G ^< 2n + *)) / dkA(k)e-^^-Q\ (3.8) 

where G_|_ arises in factoring the kernel of the integral equation, (|3.6| ) : 

, . \P2r: ,—iu + E.-iu> 

G+(u) = — f — — . 3.9 

r (|-^) 27re 1 } 

We include above an extra factor of e omitted from |28] through a typo. Q is determined 
implicitly by the equation 
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This differs from [28] by a factor of 2. This same factor of 2 is missing from eqn. 8.2.38 of 
28[| which should, we believe, read 



1 f Q 

-(e d + U/2) = / d\a(X). 

n J-oo 



(3.11) 



case ii: In the next case the location of the Fermi surface satisfies the constraint 



a<Q<r ^w~w- (3 - 12) 



In this case / j_ is computed to be 



n dVl =2-V2+ ^(r 1 - Q) - ^0(1-' ~ Qf + Odr 1 - Qf). (3.13) 



The first two terms are found in |28| although we disagree by a factor of 2 in the term of 



0{I 1 — Q) while the remaining term was computed by the authors alone. 

case iii: In the final case we are far from the symmetric point such that (U/2+ed) 3> VUT 
and Q> I' 1 = -1(U/8T - T/(2U)). We then have instead 



1 p oo J 

-4to / —T(l/2 + w)e 2 ™^-^sm(2nw)(-)- w , (3.14) 
2ix 6 i A J w e 

with Q in this case determined by 



Q = ?* + ^log(27re<f), 

_ e d + U/2 ^ 

q* — 



V2UT 
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Figure 4: A plot of the linear response conductance at zero temperature in zero 
magnetic field. The parameters used are U = .75 and V = U/12. The dashed line 
marks out the conductivity derived from a numerical solution of rid while the solid 
line represents the closed form solution described in this section. 

In Figure 4 is plotted the linear response conductance as a function of the dot chemical 
potential, Ed > —U/2 (for Ed < —U/2 particle-hole symmetry tells us the plot is a mirror 
image about the Ed = —U/2 axis), according to this closed form solution. For the purposes 
of comparison, we also present the conductance derived from a numerical evaluation of the 
equations determining rid- The vertical lines divide the plot according to the three cases 
of the closed form solution. We see that this solution best matches the numerical solution 
in cases i and iii. We also see that the solution makes a discontinuous transition from case 
i to case ii, a consequence of the approximate nature of the solution in case ii. 

As expected the linear response conductance rises smoothly from zero at large, positive 
values of Ed to its maximum possible value, 2e 2 //i, at the symmetric point of the model, 
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U/2 = —Ed- The ratio of the values of U and V chosen for this plot correspond to that of 
the experimental realization of a quantum dot discussed in |iJ0. 

3.2. H 7^ Linear Response Conductance 



2 



1.5 



J 9 1 

o 
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e d /r 

Figure 5: A plot of the linear response conductance at zero temperature for various 
values of the magnetic field. By particle-hole symmetry, the conductance for values 
of Ed < —U/2 is obtained by taking the plot's mirror image about the axis Ed = 
—U/2. The parameters used are U = .75/nD (D being the bandwidth) and V = 
XJ/12. For these parameters, the Kondo temperature at the symmetric point is 
T fc = .02508r. 

Solution Away from the Symmetric Point (Ed > —U/2): 

Again the transmission amplitude of the electrons is given by ( |3.3|) and ( |3.4j) . But in 
1 Note that our definition of F is related to that of [0] by F = rjjJ/2. 
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this case 

pB rQ 

n dT = / dkp imp (k) + / dA<7 imp (A) 

J ~° JQ (3.16) 

rQ 

n d i = / dXa imp (X), 
JQ 

where pi mp and o"i mp are given by ( |2.19| ). 

In general the equations for p irap and a imv cannot be solved analytically. Therefore 
we resort to numerical solutions. In Figure 5 we plot the result. Presented there is the 
linear response conductance as a function of Ed for a variety of magnetics fields ranging 
from H = to H = 3T. 

As H is increased from zero, we see two effects: the value of Ed marking the conduc- 
tance peak shifts away from the symmetric point, Ed = —U/2{= — 61?) while the magnitude 
of the peak decreases. This is as expected. The Kondo temperature for the model is given 

by m m 

Tk = y^ e -( £d ( £d +t/)-r 2 )/(2rt/) ? (3 _ 17) 

and so varies strongly as a function of the dot chemical potential. When H > Tk we expect 
the Kondo effect to be suppressed and any consequent enhancement in G to disappear. For 
values of Ed away from the symmetric point, Tk is relatively large and thus strong fields are 
needed to suppress the conductance. Closer to the symmetric point, Tk is exponentially 
suppressed and weak fields are sufficient to destroy the Kondo effect. 

When H = 0, a conductance maximum of 2e 2 /h occurs at the symmetric point. At 
the symmetric point, fi^yi = 1/2, and so each spin species makes a corresponding e 2 /h 
contribution to G. As H is increased to large positive values, the gate voltage, Ed, at which 
ndi /i = splits leading to a corresponding split in the conductance resonance. For 
example for large H the resonance associated with the spin j electrons is approximately 
e 2 /h and occurs at H/2. 

We see for example in Figure 5 that when H = T/lOO the Kondo temperature is 
never exceeded regardless of the value of the gate voltage and so we see little consequent 
suppression of the conductance. However for the next largest value of H, H = T/15, the 
Kondo temperature is exceeded in the Kondo regime and we see a corresponding depression 
in the conductance in this regime. For the largest value of H = 3r, we see as expected 
that the peak value is approximately e 2 /h and that it occurs roughly at Ed = H/2. 
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We note that the linear response conductance curves are symmetric about their peak 
value. This differs from the prediction based upon a Hartree-Fock computation of Ng and 
Lee . But it is in agreement with Meir and Wingreen ||. 

The conclusions in the above discussion are reiterated in Figure 6. There we plot the 
behaviour of the conductance peak as the magnetic field is increased from zero. In the top 
panel of Figure 6 we see that the location of the peak rapidly moves away from Ed = —U/2 
towards the large field value of H/2. The straight line in this panel indicates the behaviour 
of the peak if interactions were absent. Similarly we see the peak height in the middle 
panel of Figure 6 change from its maximal value of 2e 2 /h at H = to e 2 /h at large fields 
corresponding to a contribution to the conductance of a single spin species. And finally in 
the bottommost panel of Figure 6 we examine the width of the peak. At H = the width 
of the peak is approximately 12r. However in the large field limit this settles down to 2r 
appropriate to the conductance being governed by the Breit-Wigner formula 

e 2 r 2 

G= ^T 2 + (e d -H/2r (3 ' 18) 
appropriate to a single non-interacting electron species. 

Solution at the Symmetric Point, (ea = —U/2): 

Although we cannot in general express the magneto-conductance in closed form, we 
can do so at the symmetric point. At the symmetric point, the Fermi surface of the spin 
excitations, Q, goes to oo. The density equations then simplify to 



1 , f°° 
Pbuik(fc) = — + g (k) / dXa 1 (g(k) - A)cr bu ik(A); 

27T J-oo 



o"buik(A) 

and 



x'(X) 



TV 



/oo pB 
a 2 (A' - A)a bu ik(A') - / dk a x (X - g(k))p huik (k), 
-oo J-D 



(3.19) 



/oo 
d\a 1 (g(k) - A)a imp (A); 
-oo 

/oo pB 
d\a 2 (X' - A)cr imp (A / ) - / dk ai(A - g(k))p imp (k). 
-oo J-D 



(3.20) 
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Figure 6: Plots of how the conductance peak evolves with increasing magnetic field. 
In the top panel is a plot of the location of the peak while the middle panel records 
the peak height and the bottom panel gives the peak width. The parameters used 
are U — .75/7rD (D being the bandwidth) and V = U/12. 
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Here the limit B is determined by 




(3.21) 



As the electrons in the leads are non-interacting, the first equality is a result of Pauli- 
paramagnetism. 



closed form 
numerical 




0.4 0.6 

H/T k 

Figure 7: A plot of the linear response conductance at zero temperature at the 
symmetric point {ea = —U/2) as a function of magnetic field. The parameters used 
are U = .75 and T = U/12. The dashed line marks out the conductivity derived from 
a numerical solution of while the solid line represents the closed form solution 
described in this section. 



The phase shifts are given by 



rQ rB 
<5 eT = 2% - 8 e ± = 2% d\a imp (X) + 2tv dkp imp (k) 

JQ J-oo 



71 + Ti dkp imp (k), 

J — oo 



(3.22) 
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which follows as 

.B 



/ d\a imp (X) = - - - / dkp irnp (k). (3.23) 

J -oo " " J — oo 

This is established by integrating, f_ , ( |3.20| ). We can thus focus solely upon the k- 
distribution. 

In order to evaluate the phase shift in the equation above, we are thus interested in 
computing the integral, 

/ dkp imp (k) = 2M l5 (3.24) 

J — oo 

which, as indicated, is directly related to the impurity magnetization, Mj. Using the same 
Wiener-Hopf technique, |23] determined this integral in the case Tj~ > H, to be 



2M- = — V G +( i < 2n+1 )) ( 1)n -n(2n + l)( b -l/I). 

7T ^ (2n + l) K ' 



(2n + l) v ~ J " 

(3.25) 



7 l w 2 UT 

b= - log — \ ). 

7T H V Tre 



Combining this with the expression for the Kondo temperature in (|3.17|) , we have for the 
scattering phases at leading order in H/Tk 

5 eT =27r-5 ei =7r(l + -^), (3.26) 

which in turn gives the magneto-conductance as 

G(H) = 2^(1 - ^(^) 2 + 0(^n (3.27) 
h 16 T k T k 

The quadratic deviation from the maximal conductance has the expected Fermi liquid 
form. 

In Figure 7 is plotted how the magnitude of the linear response conductance at the 
symmetric point chang function of H/Tk (for small H) according to this solution. 

We plot it against the numerical solution and we obtain agreement at worst of 1.5%. The 
disagreement becomes larger as H increases as the closed form solution is only valid at 
H < Tk oc e -1 /^. Beyond H > T^, we rely entirely on a numerical solution to determine 
the magneto-conductance. We plot the behaviour of G(H) in Figure 8 up to H = T. 
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Figure 8: A plot of the linear response conductance at zero temperature at the 
symmetric point (ed = —U/2) as a function of magnetic field. The parameters used 
are U = .75 and V = U/12. The dashed line marks out the conductivity derived from 
a numerical solution of rid while the solid line represents the closed form solution 
described in this section. 



4. Linear Response Conductance at Finite Temperature 

In the previous section we focused upon scattering at the Fermi surface. In this 
section, we discuss a problem that requires us to understand scattering at finite energy: 
the linear response conductance as a function of temperature. We compare it to Costi 
et al.'s NRG results and find excellent agreement. This is important as it indicates we 
have an essentially correct description of the low energy scattering states, i.e. an excellent 
approximation to the right hand side of ( 2.52Q . 

Computing the linear response conductance at finite T is a complicated matter. Even 
though we only consider glued charge and spin excitations as explained in Section 2, we now 
have to compute their scattering matrix - the 1/L correction of their associated momenta - 
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in the presence of a "thermalized ground state" . This ground state is no longer composed 
of merely real k states and 2-string bound states of spin and charge as it was at T = 0. 
Rather all the possible solutions of the Bethe ansatz equations of the model make an 
appearance. Thus to begin the computation of the scattering amplitudes at finite T, we 
give the complete list of excitations in the model. 

It is useful to understand how the following calculation differs from the exact com- 



putation of the conductances in the fractional quantum Hall problem ||24|| . The logic of 
2^] would first require the identification of all the excitations above the zero temperature 



ground state. We then would need to compute both the bulk and impurity scattering 
matrices of these excitations. Having done this, the second step would be to turn on the 
temperature. The thermalized ground state would them consist of seas of all these exci- 
tations with their distributions determined by a thermodynamic Bethe ansatz, employing 
the zero temperature S-matrices. Finally, the conductance would be determined by the 



zero temperature impurity scattering matrices as in ||24fl . A potential difficulty in this 
approach would come in understanding the gluing necessary to form excitations carrying 
electronic quantum numbers as it only such excitations that can be mapped back to the 
two-lead picture. 

Whatever the likelihood for success, one can see that this way of proceeding requires 
considerable technical expenditure. It turns out to be easier to determine the thermal- 
ized ground state directly from the Bethe ansatz, i.e., by dealing with bare excitations 
rather than first determining the zero temperature ground state, classifying its "physical" 
or "dressed" excitations, and expressing the thermalized ground state in terms of these 
excitations. On the other hand, with the quantum Hall edges, it was the natural way to 
proceed as the initial data one had at hand (with no calculations whatsoever) is precisely 
the excitations about the zero temperature ground state together with their S-matrices, 
both bulk and impurity. 

To proceed then, the first step is to identify all solutions of the Bethe ansatz equa- 
tions. These are of three types: real k-states, spin complexes associated with complex k's 
(both of which we have seen as the ground state at zero temperature is composed of such 
excitations) and spin complexes not associated with complex k's. Below we give a more 
specific description together with the quantum numbers carried by each excitation. The 
spin quantum numbers are measured relative to a vacuum carrying spin N/2 where N is 
the number of particles in the system. 
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i) real k's: These appear in the ground state at T=0 in the presence of a magnetic field. 
They carry charge e and no spin. 



ii) n-spin complex with no associated k's: An n-complex involves n X's organized as 



77,-1-1 

X n3 =X n + i(^-j), j = l,---,n. (4.1) 



Here A n is a real rapidity and is known as the centre of the complex. The n-spin complex 
carries spin —n/2. 

iii) n-spin complex with 2n associated complex k's: The n-complex is organized as before 



77-1-1 

A^=A n + ;(^--i), i = l,...,n. (4.2) 



but now there are two k's associated with each A nj 



71 

g(k+^) = X n +i(- + l-j) 



71 

g(k- n3 ) = x n + i(--j). 



3 = !,-••, n; (4-3) 



These excitations carry charge 2ne and spin —n/2. The simplest of these excitations (n=l) 
appear in the ground state at zero temperature. Now the claim is that these are solutions 
to the Bethe ansatz equations ( |2.7| ) and (|2.8|) and indeed they are in the thermodynamic 



limit (the object of our concern) [^8[ . 

We can derive equations constraining the particle/hole densities of these various ex- 
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citations in the same fashion we arrive at ( p.ll|) from ( |2.7|) and (2^8). The result is 



p p (k) +p h (k) 



r + / d\a n (g(k) - \)(a pn (\) + a' pn (\)); 

^ n=l J -°° 



oo 



VhnW = + ~ I a n (X - g(k))p p (k) 



3G 



00 poo 

/ d\'A nm (\ - X')a pm (X'); 

/OO °° /"OO 

a n (A - ^(fc))p p (fc) - / ^nm(A - A> pm (A'); 
-00 m=1 J-oo 



x n (A) = V2UT Re (A + z^) 1 / 2 + + e d ) 



A n (A) = --d x S n (X) = ——9a Re 5(-V2C7f (A + i^) 1 ' 2 + U/2 + e d ) 

TV TV 2 

~ 7T d * E S(-y/2Uf (A + Un - 2k)) 1 / 2 + U/2 + e d ) 
Ztv * — ' I z 
1 — 1 ^ 



+ 5(V2UT (A + -{n - 2k)) 1 / 2 + U/2 + e d )) 



The kernels in the density equations are given by 
2n 1 



a n (X) 



(4.4) 



7T (n 2 + 4A 2 )' 

min(n,m) — 1 (4.5) 

A nm (X) = S nm 5(X) + a\ n - m \ (A) + 2 a\ n - m \ + 2k(X) + a n + TO (A). 

fc=i 

Here Pp/^ is as before while a p /hn denotes the particle/hole densities of n-strings associated 
with complex k's (in Section 2 we denoted cr p /hi by cr p / h ) and &' p / hn denotes particle/hole 
densities of n-strings not so associated. 

As stated in the introduction of this section, we construct the electronic excitations 
in the same fashion as at zero temperature, the only difference being that the excitations 
are now over the thermal ground state, not the T = ground state. In order to describe 
the scattering amplitudes we thus need to specify the impurity momentum of the p{k) 
and o"i(A) excitations. The finite temperature momenta for such excitations are as follows 
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(compare (|2.20| )): 



00 /"OO 

p^(k) = 5(h) + J2 / (Ug(k) - A) - 2tt)«p(A) + ^p(A)); 

/OO 
dkp™v(k)(6 1 (\-g(k))-2n) (4.6) 
-00 

00 „ 

+ rfA'(E lm (A-A')KT(A / ), 

m=l 

with E nm given by 

min(n,m) — 1 

S nm (A) = (0|„_ m |(A) -2n) + 2 ]T (#| n _ m | +2fc (A) - 2tt) + (# n+m (A) - 2tt). (4.7) 

fc=i 

We again can read off the contributions to the momenta and the densities and arrive at 
the all important relations: 

a tf imp (fc) = 2V mp (fc); 

(4-8) 

a A pi mP (A) = -27ra im P(A), 

still valid at finite temperatures. The scattering phases, as they are given (as in Section 
2) from the impurity momenta, p imp (k) and p 1 ] mp (A), can be computed from a knowledge 
of p imp and <j^ mp . For example, a spin up excitation constructed from a charge excitation, 
k, and a n = 1 spin-charge complex, A, has a scattering phase of the form 



= 2tt f dkp imp (k) + 2tt [ Q dA'o-i mp (A') 

J-D J\ 



(4.9) 



identical to (|2.23| ). 



Another piece of the prescription of computing the scattering amplitudes at finite 
temperature is the energy associated with the charge and spin excitations. The energies 
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of the various excitations can be derived as in Section 2 with the result 

/(-4(A)). 



f f(—f 
e(k) = k + Tj2 / d\log(^— 

n=1 J-oo J[ t 



)a n (X - g(k))] 



-£n(A)) 

Iog(/(e„(A))) = -=*„(A) - / <%'(&) log(f(-e(k)))a n (g(k) - A) 

J —oo 

oo „ 

+ ]T / dA'A nm (A - A') log(/(-e m (A'))); (4.10) 

m=l ^ 

/oo 
^'(^(/(-e^KO^-A) 
-OO 
OO „ 

+ E / d\'A nm {\-\')\og{f{-e' m {\'))), 

where f(e) = (1 + exp(e/T)) _1 is the Fermi distribution. These equations are arrived at 
by relating the densities to the energies via 

exp(e(k)/T) = p h (k)/p p (k); 
exp(£ n (A)/T) = a nh (X)/a np (X); (4.11) 
exp«(A)/T) = a' nh (X)/a' np (X). 

This relation is chosen so that the energies determine the particle/hole distributions in the 
same fashion that they do in the case of non-interacting fermionic particles, i.e., 

p p (k) = (p p (k)+p h (k))f(e(k)), (4.12) 

and likewise for u np /h and ^' np ^ h . This definition is completely consistent with that at 
zero temperature. Taking T — ► in the above recovers (|2.36|) . This is a general feature of 



energy functionals in a thermodynamic Bethe ansatz analysis. However here the energies 
are related to the densities in an additional way indicative that the bulk of the system (i.e. 
the leads) is indeed non-interacting: 

P P /h(k) = —d k e{k)f{±s(k)); 
° np/h {X) = ~d x e n (\)f(±e n (\)); (4.13) 

Z7T 

< p/h W = 7^d x e' n (X)f(±e' n (X)). 

Having specified the energy functionals, we can now determine the particular k and A 
we choose in creating an electron. As H = 0, we will simply focus on spin j electrons. In 
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forming a spin j electron we add a k-excitation in p(k) and a A hole in ai(A). The energy 
of the electron is then 

e el = e(k)-e(X). (4.14) 

We again will only allow k to vary in varying e e \ while fixing A to some A D . While at 
T = we fixed A to be Q, its value at the Fermi surface, this is not appropriate at finite 
temperature as the Fermi surface has become blurred. However we have another way to 
characterize the correct choice for A, at least at the symmetric point, which we give in the 
following subsection. 

We are now ready to specify the final equation governing the finite temperature lin- 
ear response conductance. Given that we construct the electronic excitations by gluing 
together a fixed spin excitation, A D , and a range of k excitations, these excitations are 
distributed according to the Fermi distribution, as they must be. Thus the conductance 
at finite T is given by 

,'2 



G(T)=2^ / de el (-d £el f(e el ))\T(e el )\ 2 ; 



(4.15) 

\T(e el )\ 2 = sm 2 (^5 el (e el = e(k) - e(A ), T)). 

Here, the first formula is the standard Landauer-Biittiker formula applied to the electronic 
excitations discussed in Section 2, while the second formula follows from the expression of 
\T(e e i)\ 2 in terms of phase shifts in the even and odd leads, and S e i is given by ( f4.9|) with 
A = A D . Finally, we have used the key result that the density of states (per unit of energy) 
for the electronic excitations is a constant as follows from ( f4.13| ). 

4-1. Computation at the Symmetric Point 

So far we have discussed the computation of G{T) in general terms. In this section 
we specialize to the symmetric point {ea = —{7/2). There are two reasons to do so. 
At the symmetric point the equations become more amenable to analysis. However more 
importantly, it is only at the symmetric point that we are able to compute the conductance, 
for it is only at this point that we can compute electron scattering for arbitrary energy, as 
required by ( |4.15[ ). 

The problem is a technical one. The energy functional, e(k) is bounded below. For 
example, as the analysis of |3^] []33| and |54| shows, at the symmetric point e(k) satisfies 



e{k) > -Tlog(3). (4.16) 
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Thus we are unable to compute directly electron scattering phases for energies below 
e(k = —D) — e(X = X a ). Rather for energies below this, we must compute hole scattering 
and relate this to particle scattering via 

Sl(e e i, e d ) = Sl Q (e ho = e eh -U - e d ), (4.17) 

valid when H = 0. In order to exploit then this relation, we need Ed = —U — Ed, i.e. 
Ed = —U/2. To compute <5^ D , we remove a k and a A-hole. Thus akin to (|2.45|) we have 



Sl (E h0 ) = 2n [ dk'p imp (k') + 2n f dX' a imp (X') , e(k) = -e ho + e(A = A D ), (4.18) 

J-D J\ 

with the difference, X ^ Q and e(X ) ^ 0. 

With this, we can now consider the simplifications in the structure of the equations 
that arise at the symmetric point. Following |[28| , we can recast the equations defining the 
energy functionals in a universal form for energies comparable to T^, the Kondo tempera- 
ture. Define 

1 1 2A 

0„(A) = -£n(A - - log( — )); 

1 TV 1 

1 124 

(p'Mk)) = --E(-g(k) + - \og( — )); 

1 1 2A ^ 

foW = ™<(-A + -iogb )); 

1 7T 1 



V2UT 
2n 

In our definition of 0^, we rely upon the fact that in the energy range we are interested 
in, £(k) depends solely upon g(k). With these definitions, 4> n and 4>' n satisfy the equations 

/oo 
d\'s(X - A') log(/(re„-i(A))/(n„ +1 (A))) - 5 nl e*\ (4.20) 
-oo 

where £ n = <j) n or <j)' n and s(X) = cosh _1 (7rA)/2. These equations have been analyzed by 
|32| [^3| and [^4 1 . In practice they are highly accurate in determining energies up to scales 



of 10's of Tfc's. This is a consequence of two scales existing in the problem, and y/UT. 
These equations focus on the first scale while throwing out information on the second. But 
because Tk <^ y/UT in the Kondo regime, this approximation is extremely good. 

With these equations in hand, we can determine the choice of A D . At the symmetric 
point we expect the scattering phase to be symmetric in energy, i.e. 

8 el (e)=8 el (-e), (4.21) 
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(4.22) 



regardless of the temperature. We thus fix A = A D such that ( |4.21| ) is satisfied. 

We now derive the specific equations for the impurity densities at the symmetric point. 
These equations have the initial form 

00 poo 

p^(k)=A(k)+g'(k)J2 / d\a n (g(k) - A)« p + ^)(A); 

n=l J ~°° 

/oo 
dkp™P(k)a n (\-g(k)) 
-00 

/oo 00 
dA^A nm (A-A>; m m p (A'); 
-°° m=l 

/OO POO 00 

dkp^(k)a n (X - g(k)) - d\'J2 A nm{X " 
-00 J-00 m=1 

We will recast these equations in a simpler form. We use the inverse of the matrix A nm , 

A~lW = <W(A) - s(X)(S nm+1 + «y nm _i), (4.23) 
together with the equalities 

/oo 
dX'A-^(X~X')a m (X' -X"), 
-00 

/oo 
dkA(k)a n (X-g(k)), 
-00 

to rewrite Q4.22|) as 

/oo 
dAs(A-^(fc))Ai(A) 
-00 

/oo 
dA S (A-^))K7(A) + a;7 P (A)); 
-00 



(4.24) 



/oo 
dA'*(A- A') (*ZlxW + (AO) 
-00 

/oo 
dfcp^AOaCA-^fc)); 
-00 

/oo 
dA' S (A - A') {alZi M + °t^i (*')) 
-00 



(4.25) 



-00 



p 

00 



+ S nl I dkp™P(k)s(X-g(k)). 
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We can further simplify these equations. For energies <C y/UT, it is an excellent approxi- 
mation to take 

-^ coshK^)-^)) - ^ + C ~ '<*»*^ (4 26) 

8r 2LT" 

Together with this approximation we can take 

<T = o; 

(4.27) 

<C P = 0, m > 1. 



These densities are identically zero at zero temperature and are governed by the energy 
scale y/UT. Since we work at temperatures far below this scale, they can be safely approx- 
imated as zero. With this, the density and energy equations can be solved numerically 
through iteration. 

We do so and plot in Figure 9, G as a function of T/Tk. Comparing to the NRG 
computation of Costi et al., we find excellent agreement for energies up to several T^, 
the regime where one would expect the NRG, by its very nature, to be most robust. We 
emphasize that this agreement is achieved with no fitting parameters. Our definition of 
the Kondo temperature, T^, is the same as that used by Costi et al. Because of the Fermi 
liquid nature of this problem, we know the functional form of the conductance at T <C Tk 
is 

2e 2 T 2 

G(T/T k ) = —(l-c^ + ■■■). (4.28) 
Costi et al. ||21]| . based upon results borrowed from |35|] and |36[], computed c to be 



c = — = 6.088.... (4.29) 
16 

We find numerically 

c = 6.05±.l. (4.30) 

We have arrived at this value by fitting the plot in the region T/Tk < -1- The error is 
systematic in nature, arising from the arbitrary nature of deciding the region over which 
to fit. 
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Figure 9: A plot of the scaling curve for the conductance as a function of TjT\~. 
Tk is as defined in Costi et al.'s work and so there are no free parameters. Our 
computation was carried out at the symmetric point in the Anderson model, U + 
2e d = 0. 



We also compare our results in Figure 9 with 13]. It appears that the logarithmic 
dependence in G, 

G~l/log 2 (T/T fc ), 

characteristic of weak coupling and arising from a one- loop RG ||13fl . should only be ex- 
pected to become qualitatively descriptive for values of T/Tk in excess of about 20. This 
observation will play a role in our determination of the validity of our computation of the 
zero field differential conductance in the next section. 

The quality of the fit is a good indicator of the validity of our approach in the Kondo 
regime. We expect from arguments given in Section 2 that our methodology should be 
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characterized by errors of order 0(Tk/yUT) << 1 and as such we should see an exact 
match between our scaling curve and the NRG results. We are thus uncertain whether the 
slight discrepancy between ours and Costi et al.'s results at large T is a consequence of 
the some unguessed shortcoming in our approach, some problem with the NRG, or some 
difficulty with our numerics. While we cannot speak to the first two, we do note that our 
handling of the numerics opens up the possibility for error at large T/Tk] the numerics are 
fashioned so to more readily reproduce the low temperature behaviour. 
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10" 3 10" 2 10"' 10° 10 1 10 2 

Figure 10: A comparison of the data from |jj with the computed scaling curve for 
the conductance. We now plot on the abscissa the ratio of the conductance with the 
maximal possible conductance. For the experimental realization in Jl| , the dot- lead 
couplings, Vi } 2, are asymmetric, and the conductance does not achieve its unitary 
maximum, 2e 2 /h. However the scaling behaviour of G/G max is expected to be the 
same. 



53 



We end this section by comparing in Figure 10 our scaling curve with the experimental 
results of jl]. We see that we find excellent agreement. We point however that while we 
compute the scaling curve at the symmetric point (U/2 + ea = 0) of the Anderson model, 
the data in |l] was taken away from the symmetric point but still in the dot's Kondo regime. 
(The Kondo temperature obtains an exponentially suppressed minimum at the symmetric 
point and so is usually below the temperature that can be experimentally realized. In 
order to experimentally see Kondo physics one then must move away from the symmetric 
point.) The continuing applicability of our scaling curve suggests a certain robustness to 
the scaling behaviour. 



5. Out-of-Equilibrium Conductance 




Figure 11: A sketch of two leads attached to a quantum dot. Each lead is, as 
indicated, at differing chemical potentials. 



5.1. Basic Formalism 

In order to compute the non-equilibrium conductance we imagine placing the reservoirs 
attached to each lead at differing chemical potentials as pictured in Figure 11. This on the 
face of it poses a problem. In doing do we add a term to the Hamiltonian of the form 

H* = k J (5.D 

This term does not behave well under the map into the even/odd electron basis in as much 
as the odd electron no longer decouples. It would thus seem it is not possible to employ 
the results of the previous sections in analyzing the out-of-equilibrium system. 

However we must ask what we need to compute the non-equilibrium conductance. We 
need to know the distribution of particles in each of the two reservoirs. And we need to 
know the scattering amplitudes of said particles. For the particle distributions, we note 
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that the particles in the two reservoirs do not interact with one another. Knowledge of 
one distribution is not needed to determine the other. Thus to compute the distribution 
of particles in reservoir 1 we can set \i<i to be whatever is convenient and likewise for the 
determination of the distribution in reservoir 2. This is notably different from what occurs 
in the scattering of quasi-particles between quantum Hall edges. In the boundary sine- 
Gordon formulation of this problem, the two reservoirs, one reservoir of positive solitons 
and one of negative solitons, do interact with one another. The above device would thus 
not work in this context. 

We emphasize again that the distribution of particles we compute are not the plane 
wave modes of an electron but rather are "dressed" electrons. But they do share several 
features with plane wave electrons. Beyond carrying the same quantum numbers of elec- 
trons, they share the same constant density of states as a function of energy as the plane 
wave electron modes. Moreover their dispersion relationship is the same as the plane wave 
modes. 

Having dealt with the computation of the distributions we now turn to the scattering. 
Here there is no problem. As we are using an integrable basis, we are able to compute 
the scattering in the context of the in-equilibrium model. The scattering of the basis 
of integrable excitations is unaffected by the differing chemical potentials in the leads. 
We could, if we wished, adopt a basis of (dressed) excitations that was aware of the 
finite voltage. Although technically challenging, we could also compute the corresponding 
scattering amplitudes. Such amplitudes, however, would only differ from the original by 
an overall phase. As the conductance depends upon the absolute value of the scattering 
matrices, our answer would go unchanged. The sole consequence of \i\ 7^ ^2 is then felt in 
the distributions. In this sense then, the problem is akin to scattering between quantum 
Hall edges. 

We now proceed with the actual calculation. Here we keep \i\ constant and imagine 
varying fi 2 alone. The computation divides itself into two cases: \i\ > ^2 and ^2 > A*i- 
We examine the former first. 
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Figure 12: A sketch of the distribution of particles in the leads when fj,\ > 

The case \i\ > \i% is pictured in Figure 12. As we are at zero-temperature, particles 
will only diffuse from lead 1 to lead 2. The current is thus given by 



We have expressed the current as an integral over all energies ranging from zero to [i2 — 
We note that the above integral reflects that the density of states as a function of energy 
is constant. Our particular choice of limits in the above integral is a consequence of our 



Fermi energy of lead 1 as 0. It is only in this energy range that particles are available 
to scatter in lead 1 and that are not Pauli blocked in lead 2. We choose energies to 
parameterize the particles as opposed to either of the parameters k or A. These latter 
parameters are not convenient in determining Pauli blocking in that with \i\ ^ ^2 the 
energy functionals for the leads are also not equal, i.e. ei(k/X) 7^ S2(k/\). We note that 
even though \i\ and ^2 are bare chemical potentials, they are the correct ones to use in 
determining the current. They are not renormalized by interactions, understandably, as 
there are no interactions in the bulk.i 

2 Technically this may be seen as follows. Changing the chemical potential in a lead by A/i 
yields a change in the number of particles, AN = — . The density of states per unit energy of 
the A-particles filling the ground state is cr(e) = As each A-particle is a bound state of two 
particles, the shift in the Fermi energy induced by the change in chemical potential is precisely 




(5.2) 



conventions: our equations governing the energy of the excitations 




A/i. 
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It is worthwhile commenting on the dependence of the current J upon \x\ and H2- 
Although the limits governing the energy range of excitations contributing to transport 
are a function of the difference, \i\ — of chemical potentials, the current's dependence 
upon n\ and Hz is more complicated. This reflects the dependence of the transmission 
amplitude, T 1 ^ 2 , upon \x\ alone. In particular, in the Kondo regime of this problem, the 
Kondo temperature that appears in the expressions for the current will be a function of 
Hx — Ed and will not depend at all upon /i 2 . 

Appearing in (|5.2| ) are the electron scattering probabilities, T 1 ^ 2 , from lead 1 to lead 
2. The amplitude, T 1 ^ 2 , is given from ( |2.6|) : 



IT 1 ^ 2 ^)! 2 = sin 2 ^ 1 ^ = -e,e d - Hi))- (5.3) 

The phase for an electron below the Fermi surface, as indicated above, is computed by 
exciting the corresponding hole. As indicated in the introduction to this section, scattering 
in this case is determined solely by the dynamics in lead 1. 

To compute 5ho involves exploiting particle-hole transformations. As such it is worth- 
while to consider the cases of zero H and non-zero H separately. With H = 0, we can 
compute the scattering of a spin { hole by relating it to a spin ] electron. According to 
( [Op and ( ^4lD , we have 



8 h l o( £ ho > Hi, £d - Hi) = 8l\e e i = e ho , -U - e d + Hi) 

dkplmpik) + 2tt / dX <x\ mp (A), e 1 (k) = e ho - 
-d ' Jq 

(5.4) 

Here the energies, e e i — (ii and Eho ~ Hi, are measured relative to the Fermi surface in 

II IT 

lead 1. By the SU(2) spin symmetry we then know 8 ho {£ho) — ^h ( £ ho)- Because of the 
behaviour of Ed under a particle-hole transformation, we can only directly compute out- 
of-equilibrium conductances when Ed — Hi < ~U/2, unusual in that it is on the other side 
of the particle-hole symmetric point. 

When H is nonzero the situation is more complicated. We no longer can equate spin j 
and spin j scattering. However we now can compute spin f hole scattering directly. From 
( |2.45|) we have 



s ll( £ ho > Hi,£d ~Hi) = 27r / dk'p\ {k') + 2n [ dX (X), e l (k) = -{e ho - Hi)- 

J-D JQ 

(5.5) 
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Because the bottom bound on e 1 (k) is —H, we are limited to computing spin j hole 
scattering for energies, < e < H. For spin J, hole scattering we resort to the particle-hole 
transformation used above: 

Ki( £ho > AH 5 £ d - in) = 5e T (e e z = £ho, ~U - E d + fll) 

/k nQ 
dk'plmpik') + 2tt / dX a\ mp (A), ^(/c) = e^o - /tii. 

(5.6) 

We have no similar constraint on the energy range for spin [ scattering. But we can see 
another issue arises. We are able to compute spin f hole scattering for a dot chemical 
potential, Ed — \i\ > —U/2, while for spin j hole scattering we can only perform the 
computation for Ed — (J>i < —U/2. We are thus limited in the case of non-zero H to the 
symmetric point, Ed — fii = —U/2. But given our belief that our ansatz for the scattering 
states is only valid near the symmetric point, this constraint costs us little. 

To compute the energy functional relating the parameter k to the energy we employ 
the equations, 

e\k) = k - ^ -fix - [ Q dX E 1 (X)a 1 (X - g{k)); 
1 Jo 



e\X) = 2x(X) -2nx - / dX , E 1 {X')a 2 {X' - X) 

Jq 

+ / dk g'(k)E 1 (k)ai(g(k) - A). 

J-D 



(5.7) 



These equations are identical to those of ( |2.36|) but for the presence of The Fermi 
surfaces, Q and B, are determined as before by 



e 1 (A = Q) = 0; 
£ i {k = B)=0, 

that is, the energy functionals are defined such that the Fermi energy is always zero. 
Computing the differential conductance then amounts to computing -e<9 M2 J: 



(5.8) 



e 2 

G(/ui,/i 2 ) = -ed^Jim,^) = -j^(\ T l^ 2 ( e = Pi- (ii,fii)\ + \ T }^ 2 ( £ = P2 ~ /Ui,/"i)| 2 )- 

(5.9) 

As the particle distribution and correspondent scattering in lead 1 are, as discussed above, 
only dependent upon fj,i, G has a particularly simple form: there are no terms of the form 
d ITI 2 

U H2 I ± I ■ 
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1 2 

Figure 13: A sketch of the distribution of particles in the leads when \i\ < \i<i- 

In the second case, \i 2 > A*i (pictured in Figure 13), such terms do come into play. 
Here, the current has the form 

J(/xi,/i 2 ) = -r / de(\T^\e^ 2 )\ 2 + \Tf-\e^ 2 )\ 2 ). (5.10) 

In this case particles scatter from lead 2 to lead 1. The choice of limits in the above 
integral now reflects that the Fermi energy in lead 2 has been taken to be zero. T 2 ^ 1 
can be determined by the same set of equations, (|5.3| ) to (|5*1)|), but with the energies 
and densities defined in lead 2. The expression for the differential conductance is more 
complicated than previously as the scattering matrices are determined on the basis of 
distributions in leads 2 and so the scattering varies as \ii is varied. We thus have 

G = -ed m J= U\T^\e = lH - ^ 2 )| 2 + \T 2 ^\e = m - ^ ^ 2 )| 2 ) 

p 2 ,o (5-11) 
+ T / de{d^\T^\e,n 2 )\ 2 + d^\Tl-\e,n 2 )\ 2 ). 

Here T 2 ^ 1 is given by 

|T 2 ^| 2 = sm 2 (U ho (-e,e d - = sin 2 (^(-£, -e d + fi 2 - U)). (5.12) 

We thus see explicitly 9 Al2 |T 2_i ' 1 | 2 is non-zero. When H ^ recall we can only compute 
T-f and only at the symmetric point. Given that we are varying /i2, and so varying the 
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effective dot chemical potential, we cannot compute the differential conductance for non- 
zero H in the case, fjL 2 > Mi- Moreover we are restricted to the region where Sd,—^ < —U/2, 
also as discussed previously. 

We again comment upon the dependence of the current upon Hi and \x 2 . As with the 
case \i\ > /U2, the current is not simply a function of the difference of the two chemical 
potentials. In this case however the scattering amplitudes depend solely upon \i 2 not /xi. 
In particular in the Kondo regime, the Kondo temperature is determined by difference of 
\i 2 with the dot chemical potential Ed- 

As 8 ho is given by an equation akin to ( |5.5|) , we can compute d^ 2 (T 2 ^ 1 ! 2 to be 

d/J^ 2 ^ 1 ^, ^ 2 )| 2 = ^sin(p imp (Q) + Pim P (k))(d^ 2 p imp (Q) + d^p imp (k)); 

<9 M2 p imp (Q) = 2n d\d m (j im p(\, -e d + /z 2 - U); (5.13) 

J Q 

f k 

d^ 2 Pim P (k) = 2n dk'd^pi^k , -e d + // 2 - U), e 2 {k) 
J —00 



-£. 



From the density equations in Appendix A (see A. 6), we find that with H = 0, 9 At2 ai mp (A) 
and pi m p (fc) satisfy 

^ 2 a imp (A) =0; (5.14) 



Q 

<9 M2 p im p(/c) = d m A(k, -e d + [i 2 - U) + r-— I d\a imp (X)s(X - g(k)) 



1 

UT 

'■Q 



/V 1 r<x> 

dXa imp (X)s'(X - g(k)) - — / dk'A(k')R(g(k) - g(k')) 

/"OO 

(d k g(k)) 2 / dk'A(k')R'(g(k)-g(k')), 



(5.15) 

where here g(k) = (k + ea — ^2 + U/2) 2 /(2UT). In computing d^ 2 T, we have neglected 



contributions from d m Q. We can see from the energy equations of (|5.7|) and (|5.8| ) that 



d^Q ± 0. (5.16) 

However, Q, determining the Fermi surface relative to the bottom of the band, is reflective 
of energy scales on the order of the bandwidth whereas we consider changes in \x 2 of O(Tfc). 
Hence d^ 2 Q is negligible. 
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5.2. Differential Conductance at H = 



The differential conductance in zero field is expected to fall off rapidly with a scale, ~ 
Tfc, from its linear response value near the symmetric point of ~ 2e 2 /h. The characteristics 
of this peak are related to the peak, the Kondo resonance, in the spectral weight of the 
impurity density of states as determined by the Bethe ansatz. This is similar to the 



findings of ||[12|[10| where they cast all transport properties in terms of the impurity 



density of states though as determined by lm(dd J{ ). With the Landauer-Buttiker approach 
we have adopted, all scattering quantites are ultimately related to the equilibrium density 
of states; the non-equilibrium density of states plays not part in the computation marking 
an important difference with []T(J . At the symmetric point we are able to derive in 
the case of negative bias, fi2 > Ml) a closed form expression for the differential conductance. 
Away from the symmetric point we must rely upon a numerical solution of the associated 
integral equations. 



i) Results at the Symmetric Point 

At the symmetric point (and hence only the case \i\ > [ii assuming (i 2 is being varied), 
we derive a closed form expression for the current and differential conductance. At the 
symmetric point, p imp is given by [f28| : 



g'(k) 1 

Pimp (A: < 0) = • ' 



2 cosh(n(g(k) - I- 1 )) 

(5.17) 

-( 7 / (fc)^e-^ (fc ^ 1+2n ) / dk'e- n9 ^ 1+2 ^ReA(ik'). 

n=0 J 



In order to make use of this expression we need to parameterize k in terms of the energy 
e(k). For energies not far in excess of we find in solving ( |5.2| ) with H = 0, 



e(k <0)-m=e 1 (k<0) = ^L e -^9(k)_ ( 518 ) 

7T 
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Hence the scattering phase is given by 
SUe, e d -^ = -U/2) = <^(e, -U/2) 

/oo fk 
dXa imp (X) + 2n dk'p imp (k') 
-oo J —oo 



-oo J — oo 



2 V l + (^-m) 2 /^ 2 - 

^l + 2n y ^2UT ' J L 1 >b 

n=0 

(5.19) 



where 



f, = -Ti = lj— e *((£d-vi)(e d -m+u)-r 2 )/(2ru) 

7T 7rV 2 



The latter term in the above is negligible when (s — fJ,i) ~ Tk as Tk << yUT. 
With this we can compute the current and the differential conductance: 

J(/Xi,/i 2 ) = -2^-T fc tan~ 1 ( ^ 2 = 

e 2 1 



(5.20) 



G( Ml ,/i 2 ) = - e a M2 j(/i 2 ) = 2 

i (i + (/x 2 - vir/Tf 



The simplicity of these results is striking. In our approach, it is directly related to the 
simple form of the dressed scattering phase ( |5.19| ), which only comes about at the end of 
a complex calculation. It is not clear to us whether there is a more direct way to obtain 
the results (|5.20|) . 



We observe that no log(/x/Tfc) terms appear in the above expressions for the current 
and conductance whereas we might expect such terms for large fi/Tk . In this regime such 
terms appear in weak coupling perturbative computations |TB[. However we have already 
established with our finite temperature calculation that weak coupling perturbation theory 
is not even qualitatively accurate until one exceeds scales of T /Tk ~ 20. We expect the dif- 
ferential conductance to be governed by similar considerations. Correspondingly, we would 
conclude that our scattering ansatz as applied to the zero field differential conductance is 
at least good for energies up to T/Tk ~ 20. 

Given that we are at the symmetric point, we would expect to be able to make contact 
with low energy scattering in the Kondo model as this model should produce identical 
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results to the Anderson model in the low energy regime. At low energies we have 

5 i ( _ = _ u/2) = 3 _ in -i , l-(e-/*i) a /Jg v 

2 4 + ( e -Mi) a /2? ; (5.21) 

= 7r + 2tan- 1 ((£-^ 1 )/T fc ). 



This latter form is identical to that found for spin excitations in the Kondo model |28] . In 
the exact solution of the Kondo model, the role assigned to 'charge' and 'spin' excitations 
differs from that of the Anderson model. In the Kondo model the charge excitations are 
non-interacting and so variations in the scattering phase occur solely because of changes 
in the spin sector. In this sense it is not surprising that we find the scattering phase 
of electronic excitations in the Anderson model is equal to the scattering phase of spin 
excitations in the Kondo model. If we were to compute transport properties directly in 
a two lead Kondo model, this equivalence suggests how we would have to formulate the 
scattering ansatz that governs the gluing together of excitations from the two sectors (in 
the case of the Anderson model, this is discussed in detail in Section 2). To compute the 
finite energy scattering phase for the case of the Kondo model, we would leave k at its 
Fermi surface value while varying A, the exact opposite of what we find in the Anderson 
model. 

It is also instructive to recast the impurity density of states so that it is a function of 
energy: 

Pim P (e) = -L (5.22) 

nT kl + £ 2/ Tk 

We see then that the impurity density of states is sharply peaked about zero energy with 
a peak height proportional to 1/Tfe. The spectral density of states as determined from the 
dot correlator, lm(dd'), is also sharply peaked around zero energy. In contrast however, its 
peak height is proportional to 1/T, a wildly different energy scale than the one governing 
the Bethe ansatz impurity density of states. The two quantities then are clearly different 
thus undermining an important premise of ||17|| . 



ii) Results Away From the Symmetric Point 

Pictured below is a plot of the differential conductance in zero magnetic field. We see 
that the expected qualitative features appear, namely the differential conductance sharply 
varies on energy scales related to the Kondo temperature, Tfc. 
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Although we are not exactly at the symmetric point, we must remain close in order 
to keep true to our methodology of identifying scattering states. If, for example, we were 
to compute the differential conductance in the mixed valence regime of the Anderson 
model, we would find our results unphysical. Our construction of the scattering states was 
predicated on the knowledge that in the Kondo regime the scattering phase varies on the 
smallest scale in the problem, the Kondo temperature, and that in turn, only the impurity 
density of states for the /c-excitations is governed by this scale. In the mixed valence regime 
all of these assumptions breakdown. 




-4-3-2-10 1 2 3 4 

(n-m)/T k 

Figure 14: A plot of the differential conductance in zero magnetic field. The value 
of the parameters used in the plot are V = .05, U = 10T, and Ed = — 5. 217, where 
we set the energy scale in terms of the bandwidth, D = tt. 

We see that the differential conductance curve in Figure 14 is asymmetric about 
^2 — ^1 = 0. This is a consequence of the asymmetry introduced by variations in the 



64 



effective Kondo temperature. In the regime [i\ > Tk does not vary as it is solely a 
function of [i\ — Ed and we have assumed only \ii is changing. However in the regime 
/i2 > jUi, Tfc now depends upon /i 2 — Ed and so changes in // 2 lead to changes in T^. Hence 
the asymmetry in behaviour. 



5. 3. Results at the Symmetric Point for H ^ 



At the symmetric point we develop closed form expressions for the differential conduc- 
tance. As stated previously, the nature of our construction of the scattering states suggests 
that our results for the differential magneto-conductance to become exact as H/Tk becomes 
large. Our results are in rough accordance with |TD[ - we find that for fields, H > T^, the 



H=0 differential conductance peak at zero bias divides in two, one peak for each spin 
species. Roughly speaking, the origin of the split in the differential conductance arises 
from a similar bifurcation in the impurity density of states. The spectral weight of the 
Kondo resonance present at uj = when H = divides into two resonances near eV ~ ±H, 
again one associated with each spin species. The peak at negative bias, < A*i, corre- 
sponds to a spin j resonance while the peak occurring with \ii > \i\ is associated with a 
spin I resonance. Given our ability to work only at \i2 < Hi at the symmetric point, we 
explore the former alone. 

For fields H ^> Tk, this peak is found at a bias close to —H. However unlike |10|] we 



find the differential conductance peak does not occur exactly at eV = —H. This is not 
surprising as this result was predicated upon a second order perturbative computation. 
We find instead that the peak is shifted to values of e\V\ smaller than H. For large 
fields we can develop closed form expressions for the position, height, and width of the 
conductance peaks. A related computation was done in |L7j in the context of the Kondo 
model. However there the analysis was restricted to the peaks in the equilibrium impurity 
density of states as determined by the Bethe ansatz and not the conductance per se. As 
discussed previously, the two are not directly or obviously related, as indeed is clear from 
the work here. 

In order to proceed with the computation, we review the constituent elements. The 
scattering phase for spin j hole scattering is given by 

SL(e ho > m) = 2n [ dk'p} mp (k') + 2n ^(A), e^k) = -(e ho - m). (5.23) 

J-D J-oo 
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Using (|3.23|) , we can write the phase solely in terms of Pi mp : 

/k r-B 
dk' P L v (k') +tt(1- / dfcpLpW). ^(k) = -{e h o-^i). (5.24) 
-D J -D 

The scattering phase for spin 1 hole scattering is found to be 

/A; /.B 
dk'pl mp (k')+n(l~ dfcpLpW), e 1 (fe) = e/»o-/xi, (5.25) 

through a particle-hole transformation. 

For iif satisfying T^, we can arrive at a closed form expression for the differential 
magneto-conductance. With H -C T^, the impurity density of states for the k-excitations 
retains its zero field form, 

k 1 

p imp (k) = -Jjf 2cosHn{g{k) _ I -i ) y ( 5 - 26 ) 

The impurity density is unperturbed by the field at first approximation as its spectral 
weight is found at the scale, T^, while the presence of the field only affects energies far 
below this by assumption. On the other hand the energy is shifted by a constant from its 
zero field value (again taking pi =0): 

^— (5-27) 

IT 

For k ^> B the energy is shifted by the bare energy of a spin in a magnetic field, H/2. 
For k <C B, the effect of the field upon e 1 (k) can be determined by rewriting the energy, 
e 1 (k) — > e 1 (k) —H and substituting into (A. 11). We then find that the field, if, disappears 
from the the equation, leaving us to conclude that energy is shifted by a dressed magnetic 
energy, H. 

Using these forms for the energy and the impurity density, the spin j scattering phase 
reduces to 

8Ue h0 > p,) = 5 -n- sin" 1 ^-^-Pi-Hf/Tj 1 x A — H 2 /T% ^ 
h ° 4 h+ieno-n-H)*/!*' 2 \ + H 2 /T 2 k h 1 J 

while for spin j scattering, we have 

Kn( e ho > Ui) = -r7T - sin 1 ( — — — ; y / J° ) + - sin 1 ( /v J^ ) . 

(5.29) 
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With this we can compute the differential conductance 

1 + 



,2 



e~ . 1 1 + (# 2 - (/*a ~ Mi) 2 )/T, 



2 //t (',, ,, /^p2 



1 l+jg 2 /4-(^ 2 - m ) 2 )/T 2 

2 (1 + #2/411)1/2(1 + ( M2 _ ^ _ H/2Y/TI) 

f, - ljEL e ^((e d -^i)(ed-^i+U)-r 2 )/(2TU) 

k ~ ttV 2 



(5.30) 



Taking il — > recovers ( j5.20|) . 



For values of if > Tfc, we must resort to a Wiener- Hopf analysis of the scattering 
phases. The details have been relegated to Appendix D where exact forms of the scattering 
phase and the energy, e(k), can be found. Here however we summarize their asymptotic 
forms. For k <^ B and H 3> Tk, the energy e 1 (k) as given in (D.8) takes the form 

1 1 



£ 1 (k) = -H(l- tan 



1 1 r^(l/2) , /IN t, 1 

4^ i + ( ,( fe) -6)^ rTl72T + 1 - ( ^- 6)tan W 



+ C + i log(47r 2 (l + - 6) 2 ))]) ( 5 - 31 ) 



-° Vv^ 1 + (</(*) -&) 2 g(k)-b< 



+ 

n 



+ 0((g(k) - b)'% 
where C = .577216 ... is Euler's constant and b is given by 



1 2 TTV 
6=-log(-J — )• (5.32) 
7T H \ Tie 

Note that only the first two terms in the above expansion are at leading order but we need 
include the remaining terms in order to obtain reasonable estimates of the properties of 
the conductance peak. Under similar conditions for k and H we obtain an expression for 
J dkp im p- 

2tt / dkp imp = 7r + 2tan- 1 (2(/- 1 - g(k))); 

(5.33) 

sr 217' 

where again J -1 determines the Kondo temperature, Tk ~ e~ nI . 
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Combining this analysis with numerical work and the results in ( |5.30| ) allow us to 
plot in Figure 15 the magneto-conductance for a variety of values of H/Tk. As explained 
earlier, we are able only able to compute the magneto-conductance for biases satisfying 
Hi < fj,\. Nonetheless we expect the differential conductance to be roughly symmetric 
about the V = axis. 

We see that for H ^> Tk the behaviour of the differential conductance is in rough 



accordance with the predictions of Wingreen and Meir |10[, that is, there is a peak ap- 
proximately at H2 — P>i ~ H for H ^> Tj~. However, as we have already noted, the peak 



occurs at a bias noticeably smaller than H while flTC§ finds the peak to occur exactly at 
H. As H is decreased this peak moves to smaller ratios of (p± — P2)/H before disappear- 
ing altogether at H ~ 2V For values of H -C Tk, the differential conductance does not 
appreciably change from its linear response value for voltage biases of the same order of 
magnitude as H. We also see that as H is increased the width of the peak narrows and the 
height of the peak approaches the value of e 2 /h indicating that only a single spin species 
(in this case spin j) is contributing to the conductance. 

We now analyze the properties of the differential conductance peak for values of H 3> 
Tfc. As we have already noted in Section 2, we expect our results to be exact in the 
regime that H/T^ becomes asymptotically large. In this regime the differential magneto- 
conductance is determined solely by the spin j hole scattering which we are exactly able 
to determine: no scattering ansatz is needed here. 

For H/Tk ^ 1, we can write the scattering phases as follows: 

/k pB 
dk'pl^k') + tt(1 - / dkp\ mv (k)) 
-D J-D 

w 7r + 2tan- 1 (2(/- 1 - g(k))); 



/k pB 
dk'pL^k') + tt(1 - / dkpl mp (k)) 
-D J-D 



(5.34) 



w -7r + tan- 1 (2(/- 1 -&)). 

In the first of the above equations, the second integral has been neglected relative to the 
first, valid for H ^> T^. We also make the approximation that the spin j, conductance varies 
inappreciably from its Fermi surface value as the bias is varied. Thus we set k = B in the 
second of the equations. As the conductance of the spin { species for large eV ~ H 3> 
is constant, the peak maximum occurs when 

S ho( £ ho > Vl) = 7T- (5-35) 
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This is turn implies the condition 



g(k)=I-\ (5.36) 

and so from ( |5.31| ) the bias at which the maximum occurs is 

eV max = e 1 (k = -V2UTI- 1 ) =-Il(l-±- tan" 1 - + •••), 

iff 7 "/— (5 ' 37) 

i- 1 -6 = iiog( ; fl " 




7T ov 2T fc V 2 



where • • • indicates we have not written out all the terms arising from ( |5.31|) . The half 
maxima of the peak occur when 5^ (£^ > A*i) = 17 if) which in turn implies 

g (k) = r 1 ± 1/2. (5.38) 

Hence the peak width equals 



eAV = e\k = -VWfir 1 - 1/2)) - e\k = -V2UT (J -1 + 1/2)); 

H ( _! 1 _! 1 \ (5-39) 



= — tan — = tan — = + 

Finally we can estimate the peak height. The peak maximum will be characterized by the 
maximal conductance of the spin j electrons (e 2 /h) together with the associated conduc- 
tance of the spin j electrons. The latter will vary only slightly from its Fermi surface value 
as already discussed. Hence the height of the peak is given by 



e 2 1 
C max = -(1 + sin 2 (-^ o ( £ = e F ))) 



e 



2 



3 (J" 1 - b) 



(5.40) 



h ^2 ^4(7" 1 - b) 2 + 1 
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Figure 15: A plot of the differential conductance in a magnetic field at the symmetric 
point. The value of the parameters used in the plot are V = .05 and U = 10T. 



The results for the location and width of the differential conductance peak are similar 



to those found by Moore and Wen |l7j to characterize the location and width of peaks 
appearing in the equilibrium Bethe ansatz impurity density of states for the Kondo model. 
In the large H limit, the impurity density of states as found in the Bethe ansatz then evi- 
dently shares certain properties with the non-equilibrium spectral density of states defined 
from the dot correlator, lm{dd'). However |17| makes no prediction as to the height of the 
differential conductance peak. We in general do not expect the height of the peak in the 
Bethe ansatz impurity density of states to be related to the height of the spectral density 
arising from Im^d^). We already know that no such relationship exists at H = (see 
Section 5.2i) and there is no reason to expect it to appear at finite H. 

In Figure 16 we plot how peak characteristics evolve with increasing H. For compari- 
son, we plot both the asymptotic forms ( (|5.36| )-( [5.38| )) for the peak characteristics against 
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the exact results. We see that the location of the peak slowly approaches eV = —H with 
increasing H. This approach will be logarithmic in H as 

and b ~ — log (if). Similarly the height of the peak approaches e 2 /h but again logarithmi- 
cally in H. We note that the asymptotic forms reproduce the exact results with remarkable 
accuracy even in the region, H ~ Tf~, where the assumptions of their derivation do not 
necessarily hold. 
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Figure 16: Plots describing the evolution of the differential conductance peak with 
increasing magnetic field. In the top panel is a plot of the location of the peak 
while the middle panel records the peak height and the bottom panel gives the 
peak width. The parameters used are U = .75/ttD (D being the bandwidth) and V 
= U/12. 
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6. Concluding Remarks 

Computing transport properties in a strongly interacting system is a difficult challenge. 
In this paper we have attacked the problem by combining a Landauer-Buttiker approach 
together with data from integrability. We have thus been able to provide a description of 
the scattering states in the theory that have led to several successes. We have verified the 
Friedel sum rule and provided a quantitative description of the linear response conductance 
at T=0 both in and out of a magnetic field. Our most striking result has, however, been 
the computation of the finite temperature linear response conductance scaling curve. This 
result is predicated upon an accurate description of the scattering states away from the 
Fermi surface. As such we have also been able to compute the out-of-equilibrium current, 
again both in and out of a magnetic field. In particular, we have provided a quantitative 
description of the differential magneto-conductance. Given the nature of our construction 
of the scattering states, we believe our computation of this latter quantity to become exact 
in the large field limit. 

While our technique bears a degree of resemblance to the successful, exact treatment 



of interacting quantum Hall edges |M| , the technical complexity of the two- lead Anderson 
model has prevented us from finding a definitive solution of the problem in all regimes. 
However, we have still been able to use integrability to find extremely good approximations 
to the exact results in the cases of greatest experimental interest. The situation here is 
not altogether different from the use of form-factors in calculating correlation functions: 
although the techniques of integrability do not (yet) generically lead to closed form ex- 
pressions, they are nevertheless a breakthrough, providing excellent approximations which 
are valid from the lowest energies through cross-over regimes. These approximations are 
far different from the standard, mean-field ones, for they contain all the crucial features 
of low dimensional, strong interactions. In particular the results presented here represent 
an improvement on previous approximate methods found in the literature. The quality of 
our approach can be gauged from the excellent reproduction (with no fitting parameter) of 
Costi et al.'s NRG finite temperature linear response curve. It is unlikely other analytical 
techniques could do the same. 

Our work raises two kinds of questions. The first is whether higher order approxi- 
mations can be devised in the present problem, akin to taking higher energy intermediate 
states into account in form-factors calculations of correlation functions. More precisely, is it 
possible to develop further the description of the fermions in terms of the integrable states 
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found in ( |2.52|) . The other, more practical question, is whether similar calculations can 
be performed in other models of interest, and so obtain excellent approximations of out- 
of-equilibrium transport properties based upon the exact solution of the thermodynamics 
and the proper identification of scattering states. 
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Appendix A. Practical Computation of Density and Energy Functionals 

Here we present equations for the energy and density functionals that are more 
amenable to numerical analysis. The original density equations are given by 

1 A(k) f 
P P (k) + p h (k) = — + + g'(k) / dXax(g(k) - A)a p (A); 

<t p (\) + <r h (\) = + ^ - J d\'a 2 (X' - X)a p (X') - j dk ai (X - g{k)) Pp {k). 

(Al) 

Expressing these equations in terms of the Fourier transform of cr(A) gives us 

2 \ L „ J 2n (A.2) 

a p (u) + a h (oj) = -^1 + ^1- e~^a p (uj) - J dke^ e~^ 2 Pp {k). 
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Solving for cr p (u) and substituting into the r.h.s. of both of the above equations gives 
Pp(k) + Ph(k) = ^ + ^-~ g'(k) J dXa h (X)s(X - g(k)) 

9[ } J 2tt [ n L j 2cosh(u;/2) 



-g'(k) I dk'p p (k')R(g(k)-g(k')); 

n J n L J L 

+ J d\'a h (\')R(\' - A) - J dkp p (k)s(X- g(k)), 

(A.3) 

where R(X) and s(A) are given by 



-iujX 



2tv J 1 + eM' 

1 1 

s(A) = — — — - = — I du> 



(A.4) 

1 If, e~ wX 



2cosh(7rA) 2tv J 2cosh(w/2)' 
We can further simplify the first of the equations in ( |A.3| ) by using the relations 

g '( A ) = /" rffc ^ e -to(A-g(fc)) e -|o;|/2. 

* y , f )2 ' (A.5) 

A(A) = / dfc— A(k)e-^ x -^e-^ 2 . 

Then 

+ Phik) = (i- + ^ + | cH/ifaC*) - </(*'))(^ + ^)) 

-</(*;) J dXa h (X)s(X-g(k))-g'(k) J dk' p p {k')R{g{k) - g(k% 
a p (X) + a h (X) = J dk(^- + i->(A - g(k)) 

+ J dX'a h (X')R(X' - A) - J dkp p (k)s(X- g(k)). 

These forms of the equations are better behaved numerically as 

dXR(X) = i (A.7) 
75 



(A.6) 



while 



J dXa n (X) = 1. (A.8) 

Hence when we go to solve these equations iteratively, successive iterations grow increas- 
ingly small as (l/2) n , whereas before we would not expect convergence. 

We can now derive new equations for the energy functional. As in ( |2.34| ) we find 

5E = L J dk {e+(k)5p p (k) - e-(k)5 Ph (k)} + L J dX {e + (X)5a p (X) - e~ (X)6a h {X)} 

= L J dk (k- ^)5p v (k) + 2L J dX x(X)5a p (X). 

(A.9) 

But now from ( |A.6|) we have 



S Pp (k) + Sp h (k) = -g'(k) J dX Sa h (X)s(X - g(k)) - g'(k) J dk' 5 Pp (k')R(g(k) - g(k')); 
5a p (X) + 5a h (X) = J dX' 5a h (X')R(X' - X) - J dk 5p p (k)s(X - g(k)). 



(A.10) 



Solving for 5a p /5ph and substituting into (]A.9|) leaves us with 



e+(k) +e~(k) = fc-y-2 J dXx(X)s(X - g(k)) 

+ [ dXe + (X)s(X-g(k))- [ dk' g' {k')e~ {k')R{g{k) - g{k% 

J J (A.ll) 

e + (X) + £"(A) = 2x(X) -2 dX'R(X - X')x(X') 

+ J dX'e+(X')R(X- X')+ J dkg'(k)£-(k)s(g(k)- A). 

Appendix B. Computing Scattering via the Impurity Energy 

In this appendix we compute the scattering phase of the electronic excitations through 
examining the impurity energy. To do so we will relate the impurity energy to the impurity 
momentum and then use the already established relations in Section 2. In doing so, we will 
bring out subtleties in defining the impurity energy for the purposes of deriving scattering 
phases. 

To determine the impurity energy of the excitations, we play a game similar to that 
used previously in deriving s(k)/e(X). The total impurity energy has the form (from ( |2.10| )) 

E imp = - J dX t7 p (A)(2r(A)) - J dk Pp (k)8(k), (B.l) 
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where T(A) = Re5(x(A) + iy(A)). Hence the bulk energy is given by 



-E'bulk — E — —E iuip 



1 

L 

H 5{k) 



(B.2) 



= L J dX a p (X)2(x(X) + + L J dk p p (k)(k -j+ ^ ,. 

We can thus derive equations for the bulk energy of the excitations: 

H 5(k) f 

4 u ik( fc ) + ^buikW = k -^ + ^T~ J dX £ buik(A)ai(A - g(k)); 

4 u ik(A)+e b - ulk (A) = 2(x(A) + ^) (B.3) 
" / dX'e^ lk (X')a 2 (X' - A) + J dk g '(k)e^(k) ai (g(k) - A). 

When we add particles to the system, we want to add them relative to the Fermi surface 
as determined by the host system (in other words we add them far from the impurity) . In 
light of this, we have another constraint determining the locations, Q and B, of the two 
Fermi surfaces: 

£buik(fc = B)= e bulk (A = Q) = 0. (B.4) 

Of course, any bulk quantity will not distinguish between a Fermi surface set as above or 
a Fermi surface determined by 

£ (k = B) = e(X = Q) = 0. (B.5) 

The difference between the two amounts to 1/L corrections. However this difference is 
important if one is looking at impurity quantities. See in contrast [p8| in the context of 
the Kondo problem. 

We are now in a position to relate the impurity energy to the impurity density of 
states. Differentiating the above leads to 

^buik(*0 + ^buik(*0 = 1 + ^ - 9\k) J dX d x e^ lk (X) ai (X - g{k))- 

dxe^W + ^ b " ulk (A) = 2(x'(A) + (B.6) 
- J dX'dye^ ulk (X')a 2 (X'-X)+ J dk d k e^ lk (k) ai (g(k) - A). 

Writing 

debuik = de- —de imp , (B.7) 
77 



where £i mp is the j- contribution to the energy of the excitation, and comparing to ( |2.19| ) 
leads to the relations: 



d\e- 1Ynp (\) = -9 A p imp (A) = 27rcr imp (A); 
d k £im P (k) = ~d k p imp (k) = -27rp imp (fc). 



rQ 

£im P (A) = -2n J a imp (A); 

J-D 



(B.8) 



Hence for spin f electrons, the scattering phase is given by 

Sl = -e imp (k)-e imp (\). (B.9) 

Moreover we have, 



(B.10) 



allowing us to prove the Friedel sum rule. Note that these relationships only hold due to 
our choice in defining the Fermi surface as in ( |B.4| ). 
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Appendix C. Direct Computation of the Scattering Phase 

It is possible to provide another derivation of the Friedel sum rule that involves a 
direct computation of the scattering phase (as opposed to working through the mediating 
agent of the impurity densities) . For (purely technical) simplicity we restrict ourselves to 
the case of a vanishing magnetic field where there are no real k's in the ground state. 

As we discussed in Section 2, the computation of an electron scattering phase involves 
the phases of a /c-particle and a A- hole. As we work in the zero field limit, the /c-particle 
phase is zero and we can focus solely upon the phase of the A-hole. To this end, we consider 
the bulk density of the A Q 's, cr bulk (A). In the ground state, cr bulk (A) obeys the equation 

f°° i i 1 

crbuik(A) + / d\'a 2 {\ - A')cr bu ik(A / ) = — x (A). (C.l) 

JQ n 

Following the discussion in ||37|| , the key quantity in the following will be the shift of this 
distribution when a particle or a hole is created at rapidity A. To study this quantity, we 
go back to the discrete form of the Bethe ansatz equations, which read 

M 

2nJ a + Y, 2 ( A - - A /?) = -2Lx(\a). (C.2) 

13=1 

If a hole is made at A, the rapidities shift A a — > A« and the above equation becomes 

M 

2nJ a + (a£> - A^ 1} ) - 6 2 (a« - A) = -2Lx{\$). (C.3) 

13=1 

Setting 

a bu ik(A Q ) [a« - A Q ] = ^F(A a |A), (C.4) 
it easily follows using the equation for o"b u ik(A) in ( |C.1|) that 

F(A|A)+ / (iA / a 2 (A-A , )F(A , |A) = -— 2 (A-A), (C.5) 

JQ 2n 

where 9 2 (x) = 2tan _1 (x) — n. To proceed, it is convenient to introduce the integral oper- 
ators K and L defined by 

/•OO 

K(f(X)) ee - / dA'a 2 (A - A / )/(A / ), (C.6) 

JQ 
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and 

(T-K)(I + L)=T. (C.7) 
From ( JC.7|) , and the fact that (22(A) = ^^\^2{X), it follows that 



/■oo 

F(A|A) = / dVL(A, V). (C.i 

J A 



A similar formula with a minus sign would hold if a particle were created in lieu of a hole 
at A. One can represent L as a power series if one wishes: 

l>00 

L(A,A') = -a 2 (A-A') + / dX"a 2 (X - A // )a 2 (A // - A') + (C.9) 

Hence in particular, L(A, A') = L(A', A). 

Now consider the impurity density of states, ai mp . It obeys 

cr imp (X) + / rfA'a 2 (A - X')a imp (X') = A(A), (CIO) 

JQ 

from which it follows, 

/ dXa imp (X) = -—<j ) (Q)+ dXF(X\Q)A(X), (C.ll) 
having set (/>(A) = — 2Re 5 [x(A) + iy(X)] and where we have used that </>(oo) = 0. As 

i>oo 

n d = 2 dXa imp {X), (C.12) 
JQ 



we find, 



n d = --cp(Q) + 2 / dXA(X)F(X\Q), (C.13) 
* JQ 



the key formula of this appendix. The r.h.s. of the above equation is highly suggestive: 
the first term is proportional to the bare scattering phase of the electron while the second 
term represents the dressing of the bare phase that results from the non-trivial ground 
state of the system. 

To complete this section we now explicitly demonstrate the Friedel sum rule. To do 
so we first imagine scattering the unperturbed ground state through the impurity. The 
entire scattering phase is then 

M 

X>(A/0. (C.14) 
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Now we imagine scattering the ground state plus hole through the impurity with a resultant 
total phase 

M 

-0(A) + X>(Aj). (C.15) 

0=1 

The difference of the two defines the scattering phase of the electron: 




Q 



(C.16) 



Comparing ( |C.13|) with ( |C.16| ) shows that with A = Q we arrive at 5 e i = irnd = lim^yi. 



the Friedel sum rule in the particular case when the magnetic field vanishes. 
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Appendix D. Wiener-Hopf Analysis at the Symmetric Point 

D.l. An alternative equation for e l (k) 

We first solve the equation governing e 1 (k), the energy of excitations in lead 1 relative 
to the Fermi surface. To do so we cast it in a different form than found in ( |5.7|) . For 
simplicity we assume that \i\ is zero. However finite \i\ does not change the expression for 
e 1 provided pi <C D. 

Now e 1 (k) is the energy associated with adding or removing a k-excitation. Thus 
imagine removing a k Q < B. This induces a change in the densities p(k) and cr(A). At the 
symmetric point, the unperturbed densities have the form (see ( |A.6| )): 

p(k)= Po (k)-g'(k) f dk'p(k')R(g(k)-g(k')), 

J —oo 



a(X) = a Q (X) - I dkp(k)s(X - g(k)), 

J — oo 

while the perturbed densities, p±(k) and ai(A), due to the hole at k D , are 

p x {k) = Po (k) - h(k - k ) - g'(k) f B dk'p 1 (k')R(g(k) - g(k')), 

L J-oo 

tri(A) = a (A) - I dk Pl (k)s(X - g(k)), 

J — oo 

where L is the system size. Rewriting pi by 

px{k) -> pxik) - ^-8(k - k Q ), 

yields 

Pi(k) = Po (k) + jg\k)R(g(k ) - g(k)) - g\k) f dk' Pl (k')R(g(k) - g(k'))- 

L J-oo 

1 [ B 

ai(A) = a (A) + -s(X - g{k )) - / dk Pl (k)s(X - g(k)). 

L J-oo 



(D.l) 



(D.2) 



(D.3) 



And so changes in the density, apart from the S(k — k )/L term already scaled out, are 
governed by 

Sp(k) = L( Pl (k) - p{k)) = g\k)R{g{k ) - g[k)) - g'(k) f dk'Sp(k')R(g(k) - g(k% 

J —oo 

5a(X) ee L(ai(A) - a(A)) = s(X - g{k )) - f dk 5p(k)s(X - g(k)). 

J —oo 

(DA) 
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The energy of the excitation can be expressed in terms of 8p and 5a: 

-e^K) = -{k - y) + J dk{k - H/2)5p{k) + 2 J dXx(X)5a(X) 

= - [(k - H/2) - 2 J d\x(\)s(\ - g(k ))] (D.5) 



+ 



J dk 5p(k) [k - y - 2 J dXx(X)s(X - g{k))]. 



We see that e (k Q ) depends now only upon dp. That this form for s 1 (k ) is equivalent 
to the equations in Section 2 or Appendix A can be shown using the technology found in 
Appendix C. 

Provided B < 0, we can introduce a change of variables that simplifies matters: 

p(z) = -4Q, z = g(k), k<0. (D.6) 
At energies not far in excess of the Kondo temperature we have 



(Jfe - H/2) -2 [ d\x(X)s(X - g(k)) w ^EL e -«9(k) _ H (D.7) 

J 7T 2 

The above then simplifies to 

/oo n2 

dz'Sp(z')R(z-z'), b=—; 

s\k ) = -[^ e -»<*.> - f ] + / dz8p{z)[^-e-«* - f ], (D.8) 
V :;!//r ( e --5( fc ) - 5p(u = ztt)) + —(1 - 5p(a; = 0)), 



tv v 2 

where in the last line we have expressed e 1 (k) in terms of the Fourier transform of 5 p. It 
is to the equation for 5p that we actually apply the Wiener-Hopf analysis. 

The expression for e 1 {k) is valid provided we have removed a particle, i.e. g(k ) > b 
or k < B. If we instead add a particle at k > B or g(k ) < 6, we obtain in a similar fashion 
the following description of e 1 (k): 

poo r>2 

8p(z) = R(z - g(k )) + / dz'5p(z')R(z - z'), b = — ; 

A 2UT (D.9) 

e i(fc ) = ^?£( e -^(fc ) + Sp (u = m)) - f (1 + Sp(iv = 0)). 

7T Z 
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D.2. Review of the Wiener-Hopf Analysis 



We so review the technique as presented in [g8[ on equations of the general form 

/>oo 

/(*)= / dz'f(z')h(z-z')+g(z). (D.10) 

J A 

Writing f (z) = f(z)9(±z =F A), the Fourier transform of the above equation yields 

/+(u;) + f-(u) = f + (u)h(u;) + g(u), (D.ll) 
where Fourier transforms are defined by 

a(u) = J due luJz a(z). (D.12) 

The key step in the analysis is writing 1 — h(u) as a product of functions, G^, which are 
analytic in the upper/lower planes respectively: 

l-h(u) = G +i^ G -^y (CIS) 

We can then write (p.llj ) as 

e ~ lUjA Q$) + e ~ lUjA f~^)°~^) = g{u)G-{u)e-^ A . (D.14) 
Given e~ luA f^(u) is analytic in the upper/lower half plane, applying the operators 

/ <*»' , ? + - x v (D.15) 

27TZ J UJ — [U ± id) 



to ( |D.14| ) yields solutions for both /+ and / : 



/\ —iuj' A. 



/+M = G+M _ J ^___ 9 (, 0G -(, 0e 

p wA i r i 



(D.16) 



G (U) ZTTl J UJ — [UJ — 10) 



D.3. Determination of Sp 

Applying the above analysis to ( p.8| ) , appropriate to the case of removing a particle, 
z > b (k < B), we have 



5p + (uj) = -e^ b —^ / da/ , R(cu>)G-(u>); 
2ni J u>— [uj + id) 

^^^Hff^ ^^"OTM (D ' 17) 

6= I.o g (|^). 
7T H \ Tie 
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If uj = 0, we can continue the u/-contour into the upper half plane about the branch cut 
of G~(uj) while picking up the pole at u' = 0: 

Ma , =0) =-i + * r^^^i^y r( i +(y)e -w«w-.>. (D.is) 

With uj = in ^ we find instead 

= zrr) = e"^) + e- & ^= f°° du'^^O- (-) r(± + ^"^(M-*). 

7rV2e Jo w'-i \e/ 2 

(D.19) 

If we now instead add a particle at k > B or <?(/c ) < 6, we obtain from the Wiener- 
Hopf analysis of ( p.9|) the following equations: 



8 P (uj = 0) = -±= T d 4e 2 -'^- b) yi ^V' ■ 



KV ' / ,/ , 1 p/1 , ,.,A 



(D.20) 



2e Jo w'+s r(i+w')V e 



Z}.^. Determination of 5(k) = 2n J dkp- imv 

The impurity density of states, pi mp (/c), obeys an equation of the form 

1 



Pimp(^) = y dz p imp (z')R(z - z) + - 

Pimp(z) = - ^777v- , z = g(k), 
9 (k) 



cosh(7r(^ — / x )) 



fD.21) 



provided 5 < 0, i.e if <C \/UT. As we are interested in the scattering phase, we want to 
compute 

5(k) = 2tt / dk' Pimp (k'). (D.22) 

J -oo 

If z > 6, appropriate for when we are computing the phase of a hole, 5(k) becomes 

S(k)=2ir / ^ 16 p + (u). (D.23) 

J l-K 00 — 10 

In this case Wiener- Hopf gives the solution of p + (uj) as 

2m J uj>-{u + 1 5) 



2 cosh ^ 
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Provided we assume H > T k , or roughly equivalently, I 1 > 6, this simplifies to 
8{k) = 2 tan -1 (2 (J -1 — g(k))) + n 

1 f°° e -27rw(g(k)-b) 1 

-^J duJ u r(-+c) (D25) 

roc -2^'{I- x -b) fuj'\~ U ' 1 

x dw'- ; sin(W) — T(-+u'). 

Jo u' + uj \ej v 2 

If on the other hand we are interested in the phase of an added particle (i.e. z < b), 
we compute S(k) via 

5{k) = 2tt / dk'p imp (k') 

J — oo 

rg(k) 

= 2tt - 2tt / dzp iiap (z) 

J — oo 



-oo 

= 2n-2n / : F p~(u;). 

7 2tt -zcu + 5 

The Wiener-Hopf analysis then yields for p~(u) 



p iu>b (■ i 

e / . , i , , s /\-iu>'b 



(D.26) 



This gives the scattering phase as 

^>7T 

<5(fc) = — - sin- 1 (tanh(7r(^(fc) - J" 1 ))) 

+ ^ P 7 « r(JT^) I e J (D.28) 

x P / do;'- sin(W) — r(- + u/), 

J u/ — cj V e / 2 

where P indicates the principal value of the integral is to be taken. 
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